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ABSTRACT 

In systems undergoing starbursts the evolution of the young stellar population is ex- 
pected to drive changes in the emission line properties. This evolution is usually stud- 
ied theoretically, with a combination of evolutionary synthesis models for the spectral 
energy distribution of starbursts and photoionization calculations. In this paper we 
present a more empirical approach to this issue. We apply empirical population syn- 
thesis techniques to samples of Starburst and HII galaxies in order to measure their 
evolutionary state and correlate the results with their emission line properties. A cou- 
ple of useful tools are introduced which greatly facilitate the interpretation of the 
synthesis: (1) an evolutionary diagram, whose axis are the strengths of the young, 
intermediate age and old components of the stellar population mix, and (2) the mean 
age of stars associated with the starburst, tsB- These tools are tested with grids of 
theoretical galaxy spectra and found to work very well even when only a small num- 
ber of observed properties (absorption line equivalent widths and continuum colors) 
is used in the synthesis. 

Starburst nuclei and HII galaxies are found to lie on a well defined sequence in the 
evolutionary diagram. Using the empirically defined mean starburst age in conjunction 
with emission line data we have verified that the equivalent widths of H/3 and [OIII] 
decrease for increasing tss- The same evolutionary trend was identified for line ratios 
indicative of the gas excitation, although no clear trend was identified for metal rich 
systems. All these results are in excellent agreement with long known, but little tested, 
theoretical expectations. 

Key words: galaxies: starburst - galaxies: evolution - galaxies: stellar content - ISM: 
HII Regions 



1 INTRODUCTION 

Starburst systems in the local universe, from giant HII 
regions to Starburst nuclei, are important laboratories to 
study the evolution of massive stars and physical processes 
thought to be associated with the very early stages of galaxy 
formation. These motivations, coupled to the advances in 
modeling and observational capabilities, have led to a burst 
of activity in this field during the past decade. 

By far the most common approach to infer the physi- 
cal properties of starbursts is to compare their spectral en- 
ergy distribution (SED) with models based on evolutionary 
synthesis (eg. Mas-Hesse & Kunth 1991; Olofsson 1995; Lei- 
therer et al. 1999). This technique performs ab initio calcula- 
tions of the spectral evolution of a stellar population on the 
basis of evolutionary tracks for stars covering a wide range 
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of masses, stellar spectral libraries plus prescriptions for the 
initial mass function, star formation rate and chemical evo- 
lution. The comparison between observations and models 
may focus on particular spectral features, such as stellar 
wind lines in the UV (Robert, Leitherer & Heckman 1993), 
WR features (Cervino & Mas-Hesse 1994; Schaerer & Vacca 
1998), Balmer absorption lines (Gonzalez Delgado, Leitherer 
& Heckman 1999), lines from red super giants (Garcfa- 
Vargas et al. 1997; Mayya 1997) , or on a combination of lines 
and the multiwavelength continuum (Mas-Hesse & Kunth 
1999; Langon et al. 2001). A common difficulty faced in such 
studies is the contamination of the spectrum by an under- 
lying old stellar population, which can be significant in the 
optical-near IR range. This contamination is sometimes re- 
moved by adopting a template spectrum for the old popu- 
lation (Langon etal. 2001), or else its effects are evaluated 
from the excess flux between models and data (Mas-Hesse 
& Kunth 1999). Other uncertainties include those associ- 
ated with differences between different sets of evolutionary 
tracks, incomplete or imperfect spectral libraries; stochas- 
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tic effects and numerical issues (Cervino etal. 2001); and 
possible effects of binary stars (Mas- Hesse & Cervino 1999) . 
We refer the reader to the review by Schaerer (2001) for a 
detailed discussion. 

The massive, hot stars in young starbursts photoionize 
the surrounding gas, producing an emission line spectrum 
which can also be used to constrain the properties of the 
starburst, like its age, metallicity and star-formation rate. 
In fact, emission line diagnostics of starbursts play a central 
role in this field for the simple practical reason that emis- 
sion lines are much easier to measure than stellar features. 
A great number of studies have been devoted to develop- 
ing such diagnostics through a combination of theoretical 
SEDs from evolutionary synthesis and phototoionization cal- 
culations for the corresponding nebular conditions (Garcfa- 
Vargas, Bressan & Diaz 1995; Stasihska & Leitherer 1996; 
Chariot & Longhetti 2001; Moy, Rocca-Volmerange & Fioc 
2001; Stasihska, Schaerer & Leitherer 2001; Schaerer 2000). 

Among the many diagnostic tools developed, the equiv- 
alent width of H/3 (Whs) stands out as a powerful age in- 
dicator of starbursts. The prediction, known for more than 
20 years (Dottori 1981), is that Whq decreases as a burst 
evolves. This behavior has been so extensively confirmed by 
more elaborate calculations that Wh/3 is often used as a sub- 
stitute for an age axis in studies which investigate the evo- 
lution of starbursts (e.g., Stasihska etal. 2001). Other gen- 
eral predictions of evolutionary synthesis + photoionization 
models are that the gas excitation and that the equivalent 
width of [OIIIJA5007 should decrease as a starburst evolves 
(Coppeti, Dottori & Pastoriza 1986; Stasihska & Leitherer 
1996) although such diagnostics are more critically affected 
by metallicity effects. 

These predictions, now routinely applied to infer phys- 
ical properties of starbursts, are hard to be tested directly, 
since that requires evaluating the age of a starburst without 
resorting to emission line diagnostics. In this paper we take 
a step back in time and investigate the empirical validity of 
these long known predictions by means of a simple Empir- 
ical Population Synthesis (EPS) analysis of Starburst and 
HII galaxies. This analysis allows a quantitative assessment 
of the evolutionary state of a stellar population based only 
on observed stellar features. EPS techniques have their own 
limitations (eg, Cid Fernandes etal. 2001a), but these are of 
a different nature than the uncertainties involved in evolu- 
tionary synthesis, and thus serve as an independent test of 
predictions of evolutionary synthesis models. 

Our main goals are to: 

(i) Study the stellar population properties of a large and 
varied sample of star-forming galaxies by means of an EPS 
analysis. 

(ii) Develop and test EPS-based tools to assess the evolu- 
tionary state of star-forming galaxies in a quantitative and 
easy-to-interpret manner. 

(iii) Perform a model-independent investigation of the rela- 
tion between gaseous properties and the evolutionary state 
of starbursts. 

In section ^ we present the data sets used in this study. 
Section [] deals with points (i) and (ii) above. The EPS 
method is presented and its results for star-forming galaxies 
are discussed by means of simple empirical tools designed 
to aid the interpretation of the synthesis. We also present a 



comparative study of EPS and evolutionary synthesis meth- 
ods, which serves to test and calibrate our tools to measure 
evolution. In section ^ we address point (iii) by studying 
the evolution of emission line properties of Starburst and 
HII galaxies using an EPS-based measure of the burst age. 
Section M summarizes our main results. 



2 DATA 

This investigation requires optical spectra where both 
gaseous (emission lines) and stellar properties (continuum 
and absorption lines) can be discerned. Two data sets meet- 
ing this requirement were used in this work. 

The first set, which we denote "Sample I", comes from 
the studies of Storchi-Bergmann, Kinney & Challis (1995) 
and McQuade, Kinney & Calzetti (1995), extracted from the 
atlas of Kinney etal. (1993). This sample covers both large, 
luminous galaxies with star-forming nuclei ("Starburst nu- 
clei") and smaller, weaker systems such as HII galaxies and 
blue compact dwarves. Active galaxies were discarded, with 
the exception of NGC 6221, which is dominated by Star- 
burst activity except in X-rays (Levenson etal. 2001). We 
further limit our analysis to those objects with metallicity 
estimates by Storchi-Bergmann, Calzetti & Kinney (1994), 
which leaves a total of 41 galaxies, 18 of which are classified 
as Starburst nuclei. The spectra were collected through a 
large 10" x 20" aperture, which corresponds to a circular 
aperture of 1.3 ftf 5 kpc in radius at the median distance of 
the galaxies. 

Emission line fluxes and equivalent widths (W) for this 
sample were re-measured from the publicly available original 
spectra and found to be in good agreement with those ob- 
tained by Storchi-Bergmann et al. (1995) and McQuade et al. 
(1995). All spectra were corrected for Galactic extinction us- 
ing the reddening law of Cardelli, Clayton & Mathis (1989, 
with Rv = 3.1) and the E(B — V) values from Schlegel, 
Finkbeiner & Davis (1998) as listed in NEDQ Corrections 
for internal extinction were applied based on the Ha/H/3 ra- 
tio, whose intrinsic value was taken to be 2.86, and al lowin g 
for underlying absorption components (see Section |4.1.l[ ). 
For the stellar population analysis we have measured the 
W'b of the Call K A3933, CN A4200 and G band A4301 
with respect to a pseudo continuum defined at selected pivot 
points, located at A = 3660, 3780, 4020 and 4510 A, fol- 
lowing the methodology outlined in Cid Fernandes, Storchi- 
Bergmann & Schmitt (1998). Our values agree very well with 
those of Storchi-Bergmann etal. (1995), but there were sig- 
nificant discrepancies between our measurements and those 
published by McQuade etal. (1995). 

"Sample II" comes from the Spectrophotometric At- 
las of HII galaxies of Terlevich etal. (1991), as analysed 
by Raimann etal. (2000a,b). Most of the individual spectra 
in this atlas do not have enough signal to measure stellar 
features, which prompted Raimann et al. to average them 
in order to increase the stellar signal. Out of 185 galaxies, 
they have defined 19 groups of similar characteristics. Each 

1 The NASA/IPAC Extragalactic Database (NED) is operated 
by the Jet Propulsion Laboratory, California Institute of Tech- 
nology, under contract with the National Aeronautics and Space 
Administration. 
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group is then treated as if it corresponded to an individ- 
ual galaxy. Three of these groups are composed of Seyfert 2 
galaxies. These were kept in our analysis only to illustrate 
their systematically different properties with respect to the 
remaining groups, 10 of which are composed of HII galax- 
ies, 2 are Starburst nuclei and 4 are classified as intermediate 
Hll/Staburst systems. The typical aperture covered by these 
spectra correspond to an equivalent radius of 0.8h^ 5 kpc. 

Raimann etal. (2000a) have measured absorption line 
W's and continuum fluxes for these groups following the 
same methodology as above. Emission line properties were 
analysed by Raimann et al. (2000b). Line fluxes were initially 
measured after subtraction of a stellar population model 
which included internal reddening. The resulting Ha/H/3 
ratio was used to further correct for residual extinction to- 
wards the line emitting regions. We have adopted both stel- 
lar and nebular properties as published by these authors 
without further corrections. 

In summary, of all stellar and nebular properties com- 
piled for these 2 samples, the following will be used in the 
analysis below: (1) the W's of the Call K, CN and G-band 
absorption features; (2) continuum fluxes at 3600, 4020 and 
4510 A; (3) emission line fluxes and W's of strong opti- 
cal lines: [OIL1A3727, H/3, [OIII]A5007, Ha and [NII]A6584; 
(4) nebular abundances, as listed in Storchi-Bergmann et al. 
(1994) and Raimann etal. (2000b); (5) the "activity class", 
as reported in the afore mentioned papers. This last item is 
used to distinguish small systems like HII and blue compact 
galaxies from Starburst nuclei, which live on larger and more 
luminous galaxies, usually late type spirals. The latter kind 
of galaxies present a more complex mixture of stellar popu- 
lations and are richer in heavy elements than HII galaxies, 
as will become clear in the analysis that follows. 



3 EMPIRICAL POPULATION SYNTHESIS 
ANALYSIS 

3.1 The method: input and output quantities 

In order to provide a quantitative description of the stel- 
lar populations for galaxies in Samples I and II, we have 
used their absorption line W's and continuum colors as in- 
put to the EPS algorithm developed by Cid Fernandes et al. 
(2001a). The code decomposes a spectrum onto a base of 12 
simple stellar populations of different ages and metallicities 
(Z). This base was defined by Schmidt et al. (1991) out of a 
large sample of star clusters originally observed by Bica & 
Alloin (1986a,b). The main output of the code is the pop- 
ulation vector x, whose 12 components carry the fractional 
contributions of each base element to the observed flux at a 
given normalization wavelength Ao. This vector corresponds 
to the mean solution found from a 10 s steps likelihood- 
guided Metropolis walk through the parameter space. Since 
colors are also modeled, extinction enters as an extra param- 
eter, but this will not be directly used in our analysis. Some 
EPS studies (e.g., Bica 1988) impose that solutions follow 
well behaved paths on the age-Z plane spanned by the base, 
in order to force consistency with simple scenarios for chem- 
ical evolution. Here we follow Schmidt etal. (1991) in not 
imposing such a priori constraints in order to allow for more 
general scenarios, such as systems undergoing mergers. 



There are no major conceptual differences between this 
EPS method and that originally developed by Bica (1988) 
or its variants, which have been applied to many stellar pop- 
ulation studies (e.g., Bica, Alloin & Schmidt 1990; De Mello 
etal. 1995; Kong & Cheng 1999; Schmitt, Storchi-Bergmann 
& Cid Fernandes 1999; Raimann etal. 2000a). However, in 
this work we will explore novel ways of expressing the re- 
sults of the synthesis, which use the population vector to 
construct easy-to-interpret diagrams and indices. 

The results presented below were all obtained feeding 
the EPS code with just 5 observables: The W's of Call K, 
CN and the G-band, plus the F3660/-F4020 and F4510/-F4020 
continuum colors. The errors on these quantities were fixed 
at 0.5 A for W K and W G , 1 A for Won, and 0.05 for the 
colors. As discussed by Cid Fernandes etal. (2001a), the 
combination of observational errors, little input information 
and quasi-linear dependences within the base hinders accu- 
rate estimates of all 12 components of x, but reliable results 
are obtained grouping x components of same age. We have 
therefore employed age-grouping schemes in our analysis. 

The base spans 5 logarithmicaly spaced age bins: 10 6 , 
10 7 , 10 8 , 10 9 and 10 10 yr. Components with these ages are 
combined onto a reduced 5-D population vector, whose com- 
ponents are denoted by xq, x-?, x%, x$ and x\o respectively. 
We will also work with an even further reduced descrip- 
tion of stellar populations, in which the 10 9 and 10 10 yr 
old components are re-grouped onto xo = x$ + xio, and 
the young 10 6 and 10 7 yr components are binned onto 
xy = xq + X7. Renaming the "intermediate age" 10 s yr 
component to xi = x&, we obtain a compact, 3-D version of 
the population vector: (xy,xi,xo)- 

Normalization requires that xy + xi + xo = 1, while the 
positivity constraint implies that all components are > 0. 
Therefore, an EPS solution is confined to a triangular cut 
of a plane in the (xy , xi, xo) space, which facilitates the 
visualization of results (see Cid Fernandes et al. 2001b for 
an application of this scheme to trace the evolution of cir- 
cumnuclear starbursts in active galaxies). We note that the 
description of a galaxy spectrum in terms of x-components 
is analogous to a Principal Component Analysis (e.g., Sodre 
& Stasiriska 1999), with the difference that, by construction, 
each component has a known physical meaning. 

We concentrate our analysis of the EPS results on evo- 
lutionary effects. Furthermore, we focus on the evolution of 
recent stellar generations (< 10 s yr), associated with the 
star-forming activity in Starburst and HII galaxies. Metal- 
licity effects are discussed using the nebular oxygen abun- 
dance, which essentially reflects the metallicity of the most 
recent stellar generation (Storchi-Bergmann etal. 1994). 



3.2 EPS results and the evolutionary diagram 

The results of the EPS analysis of galaxies in Samples I 
and II are listed in Table |l| both for the (x@, X7, xs,xg,xio) 
and the Young, Intermediate and Old descriptions. The nor- 
malization wavelength is Ao = 4020 A. A very convenient 
way to present the results of the synthesis is to project the 
(xy,xi,xo) vector onto the xy-xi plane. This is done in 
Figs, ^a and b for Samples I and II respectively. Dotted lines 
in these plots mark lines of constant xo, computed from the 
xy + Xi + xo = 1 condition. Note that these are actually 
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EPS Results for Samples I and II 
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GJJM504 


35 ± 9 


21 ± 10 


25 


± 


7 


6 


± 


3 


13 


± 7 


56 


± 6 


19 


± 


6 


6.9 ± 0.6 


G.UM71 


22 ± 8 


19 ± 9 


38 


± 


7 


7 


± 


4 


14 


± 7 


11 


± 6 


21 


± 


7 


7.2 ± 0.7 


G.NGC1510 


16 ± 7 


17 ± 9 


51 


± 


8 


4 


± 


3 


12 


± 7 


33 


± 6 


16 


± 


7 


7.4 ± 0.6 


G.Cam0949-2126 


20 ± 8 


21 ± 9 


30 


± 


7 


13 


± 


5 


17 


± 6 


11 


± 6 


29 


± 


6 


7.1 ± 0.7 


G.Mrk711 


24 ± 10 


38 ± 12 


15 


± 


6 


8 


± 


4 


15 


± 7 


62 


± 7 


23 


± 


6 


6.9 ± 0.7 


G.UM140 


25 ± 8 


16 ± 9 


11 


± 


7 


i 


± 


3 


11 


± 6 


11 


± 6 


15 


± 


6 


7.2 ± 0.6 


G.NGC3089 


20 ± 7 


12 ± 7 


44 


± 


8 


5 


± 


3 


19 


± 9 


32 


± 6 


21 


± 


8 


7.3 ± 0.8 


G.Mrk710 


58 ± 11 


30 ± 12 


4 


± 


3 


3 


± 


2 


5 


± 3 


87 


± 4 


8 


± 


3 


6.4 ± 0.3 


G.UM477 


27 ± 8 


20 ± 10 


31 


± 


8 


6 


± 


3 


15 


± 7 


47 


± 7 


21 


± 


7 


7.1 ± 0.7 


G.UM103 


6 ± 5 


10 ± 6 


45 


± 


8 


15 


± 


7 


23 


± 8 


17 


± 6 


38 


± 


8 


7.6 ± 1.0 


G.NGC4507 


3 ± 2 


6 ± 4 


4 


± 


3 


26 


± 


9 


61 


± 9 


9 


± 4 


87 


± 


4 


7.1 ± 2.3 


G.NGC3281 


1 ± 1 


4 ± 2 


5 


± 


3 


39 ± 11 


51 ± 11 


5 


± 3 


90 


± 


1 


7.3 ± 2.9 



Table 1. Columns 2-6: Population vector in the (xq, xj, xg, xg, xiq) description. Columns 7, 4 and 8: xy, xj and x n r espectively. All x 
components are in percentages of the total flux at Ao = 4020 A. Column 9: mean starburst age as defined in Section 3.4 Objects starting 
with a G_ are the spectral groups of Sample II. The three last entries in the table correspond to Seyfert 2 systems. 
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Figure 1. Results of the empirical population synthesis analysis for Samples I (a) and II (b), condensed on an evolutionary diagram. 
The horizontal axis xy is the fraction of light at Ao = 4020 A due to stars in the 10 6 and 10 7 yr age bins, while the fraction xj due to 
10 8 yr stars is plotted along the vertical axis. A third perpendicular axis carries the contribution xo of older populations (> 10 9 yr). 
Dotted lines indicate lines of constant xq, as labeled. HII galaxies are plotted as triangles and Starburst nuclei as squares. Filled symbols 
correspond to galaxies with WR features. Crosses in panel (b) indicate Seyfert 2s. Note that the Starburst nucleus G_Mrk710 and the 
HII galaxy G.Toll004_296, both of which show WR features, overlap at (xy,Xi) = (0.87, 0.04) in panel (b). 



straight lines in the xy-xi plane, which appear curved be- 
cause of the logarithmic axis. 

Galaxies from both samples define a smooth sequence 
from large xy to large xi, with not much spread in xo, par- 
ticularly for Sample II. The larger spread seen in Fig. [l^, 
is partly due to the fact that the data for Sample I was 
collected through apertures typically 2.6 times larger than 
for Sample II, and thus sample a more heterogeneous mix 
of stellar populations. This interpretation is supported by 
the fact that NGC 1510, which appears in both samples, 
looks somewhat younger in Sample II (see Table ^ . Another 
source of scatter in Fig. |l|a stems from the large number of 
Starburst nuclei in Sample I. These systems, represented 
by squares in both panels, live on galaxies with a signifi- 
cant old stellar component, whose effect is to drag points 
towards the bottom-left of the plot. For instance, the point 
at (xy,xx,xo) = (0.28,0.06,0.66) in Fig. ga is NGC 5996, 
whose spectrum reveals a weak starburst immersed in an 
old population (McQuade etal. 1995; Kennicutt 1992). HII 
galaxies, on the other hand, are closer to "pure starbursts". 
In fact, they were once thought to be truly young galaxies 
undergoing their first star-formation episode (Searle & Sar- 
gent 1972). Only recently it has become clear that they too 
contain old stars (Telles & Terlevich 1997; Schulte-Ladbeck 
& Crone 1998; Raimann etal. 2000a). This explains why 
Sample II, which is essentially an HII galaxy sample, exhibits 
a more well defined sequence in Fig. |l[ with all non-AGN 
sources bracketed by the xo = and 30% contours. 



The three deviant crosses spoiling the Sample II se- 
quence in Fig. |l|b are the Seyfert 2 groups, with their pre- 
dominantly old stellar populations (Raimann etal. 2000a). 
G_UM103 has a significant "post-starburst" component, 
reminiscent of more evolved starburst + Seyfert 2 composite 
systems, while the other two groups occupy a region charac- 
teristically populated by LINERs and non-composite Seyfert 
2's (Cid Fernandes etal. 2001b). 

Since the location of a galaxy in Fig. [j] reflects the evo- 
lutionary state of its stellar population, we interpret the dis- 
tribution of objects in this diagram as an evolutionary se- 
quence, with the mean stellar age running counter-clockwise. 

There are several reasons to interpret Fig. |l| as an evo- 
lutionary sequence. First, metal absorption lines become 
deeper and galaxy colors become progressively redder as one 
moves from large xy to large xi along the sequence. In fact, 
the sequence defined by Sample II follows very closely the 
blue to red (young to old) spectral sequence in Figure 1 of 
Raimann etal. (2000a). Second, all Sample II galaxies in 
which WR features have been detected (those marked by 
filled symbols in Fig. |l|b) are located in the large xy re- 
gion of the diagram, consistent with the young burst ages 
(a few Myr) implied by the mere presence of WR stars. 
Filled symbols in Fig. |j]a mark galaxies from Sample I 
which are listed in the WR-galaxy catalog maintained by 
D. Schaerer (webast.ast.obs-mip.fr/people/scharer). Their 
more even distribution, as compared to Sample II, is due 
to the old population and aperture effects discussed above 
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and nicely illustrated by Meurer (2000). Whereas the data 
analyzed here pertains to kpc-scales, spectra used to classify 
a starburst as a WR-galaxy are usually obtained through 
much narrower slits centered on the brightest cluster, thus 
favoring the detection of young systems. Processing such 
spectra through our EPS-machinery would surely move the 
filled points in Fig. |l]a towards younger ages. (Conversely, 
one would expect that narrow slit observations of galaxies 
represented by empty symbols in the bottom right of Fig. ^]a, 
such as Mrk 357 and UGCA 410, have a good chance of re- 
vealing WR features.) Finally, galaxies located in the large 
xi zone in the top-left (such as NGC 1800 in Sample I and 
group G_NGC3089 in Sample II) have spectra typical of 
a "post-starburst" population, with pronounced high order 
Balmer absorption lines typical of A stars (Gonzalez Del- 
gado etal. 1999). It is important to remark that neither the 
presence of WR features nor Balmer absorption lines were 
used in the EPS analysis, and yet the EPS results are com- 
patible with the information apported by these observables. 



3.3 EPS analysis of theoretical galaxy spectra 

A straight-forward theoretical reason to interpret Fig. [j] as 
an evolutionary sequence is that, schematicaly, a simple 
(i.e., coeval) stellar population moves on this diagram from 
(xy,xi,xo) = (1,0,0) at age t = to (0,1,0) after some 
10 8 yr and then to (0, 0, 1) for ages > 10 9 yr. 

In order to follow this evolutionary path more closely 
we have carried out an EPS analysis of theoretical galaxy 
spectra from GISSEL96, the evolutionary synthesis code of 
Bruzual & Chariot (1993). The theoretical spectra were pro- 
cessed in exactly the same way as the real spectra of Samples 
I and II. Instantaneous burst and continuous star-formation 
models were computed for various ages between t = and 15 
Gyr, a Salpeter IMF between 0.1 and 125 M and Z = Z@. 
GISSEL96 uses stellar tracks from the Padova group and 
offers a choice of spectral libraries. We have chosen the one 
which uses the atlas of Jacoby, Hunter & Christian (1984) for 
the optical range. The absorption features necessary for our 
EPS are clearly denned with this library. Spectral resolution 
was in fact the reason we have chosen GISSEL96 over the 
Starburst99 code of Leitherer etal. (1999), which is more 
taylored to study young stellar populations but currently 
works with an optical library too coarse for EPS analysis. 



3.3.1 Burst models 

The results for an instantaneous burst are shown as a solid 
line in the (xy, xi,xo) evolutionary diagram of Fig. |2]a. La- 
bels next to selected points indicate the model age in Myr. 
As expected, evolution proceeds from xy to xi to xo, such 
that a position on this diagram can be associated with an 
age. Fig. ^a provides an alternative representation of this 
diagram, in which all 3 components are explicitly plotted in 
a face-on projection of the plane containing the (xy , xi , xo) 
vector. The idea for this projection was borrowed from simi- 
lar plots by Pelat (1997, 1998) and Moultaka & Pelat (2001). 

In principle one would expect all evolution prior to ages 
of ~ 10 9 yr to progress along the xo = contour, whereas 
in practice the GISSEL96 models oscillate between xo = 10 
and 15% for t < 2 x 10 s yr. Similarly, the value of xi starts 



to increase before t = 10 7 yr, when all stars should still be- 
long to the xy age bin. These deviations occur due to the 
limited number of observables used in the synthesis and be- 
cause these contain observational errors which broaden the 
likelihood-function of x in a non-trivial way (Cid Fernandes 
et al. 2001a). As a result, some of the true xy proportion al- 
ways spills over onto xi and xo, and so on. Overall, however, 
these figures show an excellent correspondance between em- 
pirical and evolutionary populations synthesis calculations. 

In Fig. ^c we show the behavior of the 5 age components 
x^-xio as a function of the age of the GISSEL96 models. 
The population vector evolves smoothly with age, except 
for the small kink a little short of 10 7 yr due to the sudden 
appearance of red supergiants (Chariot & Bruzual 1991). 
As expected, Xe peaks around t = 10 6 yr, xj peaks around 
t = 10 7 yr and so on, but note how X7 is less well defined 
than any other component. The EPS decomposition tends 
to represent a ~ 10 7 yr burst as a combination of x§ and xg 
instead of a strong x-j. Also, a non-negligible fraction of the 
10 6 yr component spills from xa onto X7 for t ^ 10 6 yr. A 
similar effect occurs with xg and X10 for t 10 9 yr. 

Such imprecisions in the mapping between evolutionary 
and empirical population synthesis are largely suppressed 
in the coarser, but more robust, (xy , xi, xo) description, 
as shown in Fig. ^b. The figure also shows the evolution of 
xsb = xe+X7+x$, which we hereafter treat as the "starburst 
component", representing the past ~ 10 s yr of the history 
of star formation in a galaxy. This is a more reasonable 
definition for our purposes than using only the youngest, 
ionizing population (x&), since single burst models are not 
adequate to describe kpc scale regions such as those sampled 
by the observations of Samples I and II. Instead, the inner 
kpc of star-forming galaxies, particularly Starburst nuclei, 
contains a collection of many individual associations plus a 
field population with a spread in age. The detailed studies 
by Lancon et al. (2001) and Tremonti et al. (2001) illustrate 
this point (see also Calzetti 1997; Legrand etal. 2001). Such 
systems are frequently better represented by models with 
multiple bursts or continuous star formation over ~ 10 s yr 
(Meurer 2000; Meurer 1995; Coziol, Barth & Demers 1995; 
Coziol, Doyon & Demers 2001). 

3.3.2 Burst plus an underlying old population 

Real galaxies have a mixture of stellar populations of differ- 
ent ages, and galaxies in Samples I and II are no exception. 
The ongoing star-formation which makes them classifiable 
as starburst systems is observed atop an old (> 10 9 yr) stel- 
lar substrate formed in the earlier history of the galaxy. In 
our young, intermediate and old description, the effect of 
this underlying population is to dilute the values of xy and 
xi, which represent the recent history of star-formation. As 
a result, an instantaneous burst occurring on top of an old 
background does not follow the evolutionary sequence traced 
by the solid line in Fig. |^a. 

Two quantities suffice to examine these diluting effects: 
the fraction xoid{0) of the total Lsurst + Lou luminosity at 
the start of the burst (t = 0) which is due to old stars, and 
the function l(t) — LBurst(t) / LBurst(0), which describes the 
luminosity evolution of the burst in units of its initial lumi- 
nosity. Naturally, all these quantities refer to the same wave- 
length, Aq. With these definitions, and considering that Low 
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Figure 2. Evolution of GISSEL96 models on the (iy,x/,io) EPS diagram, obtained processing the theoretical galaxies as the data. 
Dotted lines indicate lines of constant xq. (a) The solid line shows results for an instantaneous burst. Numbers indicate the model age in 
Myr (note that evolution proceeds counterclockwise). Dashed lines correspond to instantaneous bursts on top of an old population which 
at t = accounts for xoid(0) = 10% (top) and 50% (bottom) of the flux at Aq = 4020 A. (b) Models with continuous star formation. 




does not evolve significantly on the time-scales of interest 
(< 10 8 yr), it is easy to show that 



7 8 9 

log model age [yr] 

Figure 4. EPS parameters against model age for instantaneous 
burst models, (a) Mean ages, (b) population vector x in the re- 
duced Young, intermediate and old description, and (c) in the 
5 component description. The t$B curve in (a) is stopped at a 
few times 10 s yr, after which it loses meaning for instantaneous 
bursts. 



Xoid(t) = 



Xoid(0) 



XBurst{0)l(t) + XQld(0) 



(1) 



where xsurst(0) = 1 — xoid(0). We can now look at the evo- 
lution of x presented above for a pure burst as corresponding 
just to the LBurst(t) component, which at time t accounts 
for only XBurst(t) = 1 — xoid(t) of the total luminosity of 
the system. This allows us to, with the help of equation [j], 
re- normalize the evolution of (xy , xi , xo) to this new scale 
for any desired value of the contrast parameter xom(0). 

Results for xoid(0) = 10 and 50% are shown as dashed 
lines in Figs, ^a and ^a. As the burst fades, the evolution- 
ary sequence bends over towards large xo quicker for larger 
ioiii(0), i.e., for smaller initial ratios of burst to underly- 
ing old population power. Despite this effect, evolution still 
proceeds in an orderly counterclockwise fashion. 



3.3.3 Continuous star-formation models 

Figs, ^b and ^b show the evolution of (xy, xi, xo) for GIS- 
SEL96 models with continuous star-formation. Since in this 
regime there are always young stars, at any given age > the 
system looks younger than an instantaneous burst model. At 
t = 10 8 yr, for instance, (xy , xi, xo) ~ (0.6,0.3,0.1) for the 
continuous star-formation models and (0.3,0.6,0.1) for an 
instantaneous burst. For t > 10 9 yr, x converges to a re- 
gion around (0.3, 0.4, 0.3), instead of plunging towards large 
xo as for instantaneous bursts. Evolutionary sequences for 
models with time-decaying star-formation rates (say, expo- 
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Figure 3. Evolution of GISSEL96 models in the xo) EPS-space seen from a face on projection of the xy + xj + xo = 1 plane, 

(a) Instantaneous burst models. The top curve corresponds to a pure burst, while the middle and bottom curves correspond to bursts 
atop a pre-existing old population with XQidi®) = 10 ana 50% respectively, (b) Continuous star formation models. 



nentially or an "extended burst" step function), would define 
curves intermediary between those traced in Figs, ^a and b. 

Since the luminosity increases without bounds for con- 
tinuous star formation, any underlying old population is 
quickly outshone by the new stars. We therefore do not 
present dilution curves such as those computed for an in- 
stantaneous burst, since, except for the earliest ages, they 



are practically identical to the undiluted curves in Figs. |2|b 
and ^b. The evolution of the population vector for contin- 
uous star-formation models is shown in Figs. [|b and c. As 
expected, the curves are smoother, and each component lives 
longer than for an instantaneous burst (Fig. m. 
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Figure 5. As Fig. Q for models with continuous star formation. 



3.4 Mean stellar age 

The experiments above demonstrate that the evolutionary 
state of a starburst can be assessed by its location on the 
(xy , xi, xo) diagram. In order to translate this location 
into a number which quantifies the "evolutionary state" one 
could use, for instance, the angle 8 = tan -1 xi/xy, which 
increases as a burst evolves. Alternatively, we may use x 
to compute the mean age t of the stellar population. Since 
stellar populations evolve in a non-linear way, it makes more 
sense to define t from the mean log t among the populations 
represented by the base: 



logi(Ao) = y^Xi(A )log 



(2) 



The dependence of x on the normalization wavelength 
is explicitly written in this equation to emphasize that this 
definition of t is A-dependent. This happens because the Xi's 
are flux fractions at Ao, so t(Ao) is a flux-weighted mean age. 
A A-dependent age makes observational sense for the sim- 
ple reason that young stars are bluer than old stars, which 
makes t an increasing function of A. Though the value of t 
depends on Ao, the evolutionary sequence traced by this in- 
dex independs on the choice of normalization, and so it can 
be used to rank populations on different evolutionary states. 

For our 5-ages base (U = 10 6 , 10 7 , 10 8 , 10 9 and 10 10 
yr), £(Ao) becomes (in yr) 



logi(Ao) = 6x6 + 7x7 + 8x8 + 9xg + lOxio 



(3) 



Since we are primarily interested in quantifying the evo- 
lutionary stage of populations associated with the most re- 
cent star-formation in starburst systems, it is interesting to 
consider a definition of t which removes the diluting effects 
of an underlying old stellar population. This can be done 
re- normalizing xe + X7 + xs to 1, which yields the following 
definition for the mean starburst age (also in yr): 



Note that by construction 6 < logi < 10 and 6 < log tss < 
8. 

The solid lines in Fig. compare our t and Jsb EPS- 
based age indices for Ao = 4020 A with the corresponding 
age of the GISSEL96 models for an instantaneous burst. 
Despite some minor oscillations, these two indices bear a ~ 
one-to-one relation with the theoretical age. Fig. ^a presents 
these same indices but for the continuous star-formation 
models. As expected, t and Jsb evolve more slowly than for 
an instantaneous burst, but they still increase steadily with 
the model age. Since t and Jsb are entirely obtained from a 
few easily measurable quantities, this result encourages their 
use as empirical clocks for stellar populations. 

As it is clear from its very definition, due to the coarse 
age-resolution of the base, our mean age index Jsb is not 
meant to be used as a fine-graded chronometer of starbursts. 
Yet, the above experiments with theoretical spectra clearly 
show that 7sb provides a useful way to rank galaxies ac- 
cording to the age of the dominant population among the 
multiple generations of stars formed in the recent history 
of star-formation. This definition is particularly well suited 
to describe spectra which are integrated over large regions 
and hence average over many such generations, Although 
all galaxies discussed here contain populations younger than 
10 7 yr, which power their emission line spectrum, this ongo- 
ing star-formation may be less intense than in the recent past 
(10 7 -10 8 yr), such that the young generations live among an 
older, non-ionizing starburst population. In this case, one 
expects to find significant X7 and xs components, and thus 
tsB > 10 7 yr. Conversely, if the current star-formation is 
more vigorous than in the past, mean ages of less than 10 7 
yr are expected. It is in this context of starbursts extended 
over a period of up to ~ 10 8 yr that we envisage t SB and 
the (xy , xi , xo) diagram as useful tracers of evolution. 

In principle, a base with a finer age resolution, including 
elements intermediate between xq, X7 and xg, could yield 
a more detailed description of the evolution of starbursts. 
In practice, however, these elements would be well approxi- 
mated by linear combinations of the existing base elements 
unless new observables were introduced in the synthesis pro- 
cess. For this reason, we opted to perform our EPS analy- 
sis with the base and observables described in 



3.1 



whose 

pros and cons have already been fully exploited in our pre- 
vious investigations (Cid Fernandes et al. 2001a,b; Schmitt 
et al. 1999). Furthermore, as we shall soon see, this relatively 
coarse description is well suited to our present purposes. 



4 THE EVOLUTION OF EMISSION LINE 
PROPERTIES 

Emission lines in star-forming galaxies are umbilicaly linked 
to their young stellar population, whose massive, hot stars 
photoionize the surrounding gas. Also, in non-instantaneous 
starbursts the continuum carries a large contribution of stars 
borne before the current generation of ionizing stars, thus 
affecting emission line equivalent widths. In this section, we 
combine the tools to measure the evolution of starbursts de- 
veloped in Section ^ with the emission line data compiled in 
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Figure 6. Equivalent width of H/3 against the EPS-defined mean 
starburst age Tsb for Sample I. Panel (a) shows the observed 



equivalent widths (Ww'|). In (b) the contribution of old stars to 



the continuum under H/3 is removed dividing Wj^S by xss(4861). 
In panel (c) Wh/3 is further corrected for the presence of an ab- 
sorption component. Triangles, circles and squares correspond to 
gas metallicity ranges of (O/H) < 0.4, 0.4-0.6 and > 0.6 (O/H)© 
respectively. Open and filled symbols are used to distinguish HII 
galaxies from Starburst nuclei. The P r values are the probabili- 
ties of no correlation in a Spearman test, small values indicating 
significant correlations. 



Section ^ to investigate whether the emission line properties 
do indeed evolve along with the burst. 



Figure 7. As Fig. ^ but for Sample II. Crosses indicate Seyfert 
2's, which were excluded from the correlation analysis. 



are determined from completely independent measurements, 
which only highlights the significance of this result. 

The different symbols in Figs. |^, Q and all subsequent 
plots represent three gas metallicity ranges, triangles, cir- 
cles and squares corresponding to (O/H) < 0.4, 0.4 to 0.6 
and > 0.6 (O/H)© respectively. Open and filled symbols are 
used to distinguish HII galaxies from Starburst nuclei. We 
postpone a d iscussion of the effects of Z and activity class to 
Sections 4.3 and 4.4, which explore emission line properties 



more directly affected by these variables. 



4.1 The equivalent width of H/3 

As a burst ages, its most massive stars are the first to die, 
resulting in a steady decline of the ionizing photon flux 
and hence on the luminosity of recombination lines such 
as H/3. The stellar continuum Chp underneath H/3 also de- 
creases, but more slowly than Lagit), since it carries a sig- 
nificant contribution from longer-lived, non-ionizing, lower 
mass stars. As a result Wnp ~ Lh/3 /Chp decreases as the 
burst evolves, as first discussed by Dottori (1981) and con- 
firmed by evolutionary synthesis calculations. Wh/3 is there- 
fore an age indicator, and it is frequently used as such in 
studies of star-forming systems (e.g, Stasihska etal. 2001). 
However, an empirical confirmation of the Wapiti) predic- 
tion requires an independent measure of the starburst age. 

In Figs. |^a and ^ we carry out this test with galaxies 
from Samples I and II respectively, using our EPS-based 
index t$B as a clock for the starburst. The anti-correlation 
between Wnp and tss is evident for both samples, thus 
confirming that Wnp decreases with time. The probability 
P r of no correlation in a Spearman's rank test is just 7x 10~ 5 
for Sample I and 10 -4 for Sample II, indicating a very high 
statistical significance. We emphasize that Wh& and tsB 



4-1.1 Corrections to Whs 

The values of Wh0 in Figs, ^a and [j]a are the raw measure- 
ments (W H b p). At least two corrections have to be consid- 
ered, both of which increase Wup- 

(1) The continuum under H/3 carries a contribution from 
an old stellar population which dilutes Wh/3 with respect to 
the value it would have in a pure starburst. Our EPS analysis 
provides a natural way of correcting for this effect, which is a 
major source of concern in studies which use Wh/3 as an age- 
indicator (Stasihska et al. 2001 and references therein). In or- 
der to isolate the contribution of the starburst to Cup it suf- 
fices to multiply it_by a;gs(4861), the starburst component 
defined in Section 3.3.1 but renormalized to A = 4861 A. The 
corrected W H p is thus simply W H b p/x SB (4861). Figs. |b and 
^3 show the dilution-corrected evolution of Wnp- 

The effects of this correction are largest for two of the 
Seyfert 2 groups in Sample II, which move significantly 
above the sequence in Fig. 0b because of their bulge dom- 
inated optical continuum (large xo)- The effect is not so 
large for GJUM103, which, as already discussed, resembles 
a starburst + Seyfert 2 composite. The fact that the ob- 
served values of Whp in Seyfert 2's are in general smaller 
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than those in starburst systems is purely due to this dilution 
effect. As explained by Cid Fernandes etal. (2001b), Seyfert 
2's ought to have intrinsically larger Wh/3 than starbursts, 
as we obtain with our EPS-based dilution correction. 

Among the star-forming galaxies in Samples I and II, 
this correction typically increases Wh/3 by ~ 50% (in good 
agreement with the corrections inferred by Mas-Hesse & 
Kunth 1999 on similar objects), but it reachs more than 
a factor of 2 in some cases. The correction is somewhat 
smaller for Sample II, partly because it contains intrinsically 
younger systems and partly because of its smaller apertures, 
which reduces "contamination" by an extended old stellar 
population. The dilution correction improves the Wh/3-Isb 
correlation for Sample I and degrades the one for Sample 
II, while for the combined sample the statistical significance 
remains unchanged at P r = 5 x 10~ 10 . 

(2) A second correction to be considered is that due 
to the presence of an absorption component hidden under- 
neath the H/3 emission. This component is present in our 
spectral base with strengths of up to Wjjb = 8 A, achieved 
for 10 8 yr populations. We have used the population vector 
x obtained in the synthesis to compute the expected value 
of W$$, typically 3-5 A. Adding to W$ yields a cor- 

rected emission Wh/3- The combined dilution and absorption 
corrected values of Wh/3 are shown in Figs, ^c and 0c. 

The absorption correction is only significant for galaxies 
with weak H/3 emission, such as those bellow Wjjp ~ 10 A 
at the bottom right of Fig. ^a. For these systems one has to 
look at the absorption-corrected values of Wh/3 as uncertain 
by as much as a factor of 2. The correction is negligible for 
most galaxies in Sample II, which, due to its objective prism 
selection criterium, contains more strong lined objects than 
Sample I. Indeed, the mean Wnp is 55 A for Sample II but 
just 22 A for Sample I. This is also why, as a whole, Sample 
II contains a higher proportion of young starbursts. 

The absorption correction degrades the Wnp-tsB cor- 
relation for Sample I slightly. For the combined Sample I + 
II data the P r value increases from 5 x 10~ 10 to 3 x 10 -8 , 
which is still significant at the 5-cr level. We thus see that 
these '1 st order corrections' introduce very little scatter, and 
do not alter our conclusion that Wh/3 does indeed evolve 
along with the stars that make up a starburst. One can also 
look at this result the other way around, and conclude that 
the fact that the Wh/3 versus Jsb diagram behaves as ex- 
pected proves the usefulness of our EPS-based evolutionary 
index Jsb, with the advantage that it is immune to the di- 
lution and absorption effects which plague Wh/3 and other 
emission-line based age indicators. 

We note in passing that these corrections alone are 
enough to bring the values of Wh/3 within the range spanned 
by evolutionary synthesis calculations such as those by Lei- 
therer etal. (1999), whereas, as it has long been known, 
the raw observed values fall bellow such predictions (Bres- 
solin, Kennicutt & Garnett 1999 and references therein). 
Differential extinction, with line emitting regions being more 
reddened than the stellar continuum (Calzetti, Kinney & 
Storchi-Bergmann 1994) and leakage of ionizing photons out 
of the HII regions associated with the starburst are further 
examples of processes that act in the sense of reducing Wh/3- 
We thus concur with Raimann et al. (2000b) and Stasinska 
et al. (2001) in that the apparent mismatch between theo- 
retical and observed values of Wh/3 bears no fundamental 
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Figure 8. Observed and predicted evolution of Wh/3- Data points 
are as in Figs, and [?]c. Lines correspond to theoretical Wupit) 
curves with t converted onto the EPS-based Jsb scale. Solid line: 
GISSEL96 models with a Salpeter IMF between 0.1 and 125 
Mq. Dashed lines: Starburst99 models with a Salpeter IMF with 
M upp = 100 Mq (long dashes) and 30 Mq (short dashes). 

physical significance. In fact, the combined effects of differ- 
ential extinction, leakage and the uncertainties in the dilu- 
tion and absorption corrections are probably responsible for 
most of the vertical scatter in Figs, pc and pp. 



4-1-2 Comparison with models 

A quantitative comparison of the predicted evolution of 
Wh[3 with that detected in Figs. ^ and [?] demands process- 
ing model spectra through the same EPS machinery used 
to analyse the data. This is necessary to translate model 
ages (t) onto our Jsb age scale. We have used the evolution 
of the rate of ionizing photons N(H°) and the continuum 
under H/3 predicted by GISSEL96 to co mpu te Wnpit) for 
the same models synthesized in Section |3.3| , for which the 
t —>tsB conversion is shown in Figs. [|a and [|a. 

The result for the continuous star formation models is 
shown as a solid line in Fig. over-plotted onto the data 
points from Samples I and II (Figs. and ^c). The dashed 
curves in this plot are the Starburst99 predictions (Leitherer 
et al. 1999) for solar metallicity, constant star formation 
models with a Salpeter IMF up to M upp = 100 and 30 M . 
These latter curves are drawn with the t -^Jsb conversion 
obtained for GISSEL96, since, as already explained, it is not 
currently possible to do an EPS analysis with Starburst99 
due to its poor spectral resolution optical libraries. Instan- 
taneous burst models (not shown for clarity) follow roughly 
the same curves up to Jsb ~ 10 6 ' 8 yr and then plunge verti- 
cally, signaling the end of the ionizing phase of the cluster. 

The rate at which Wh/3 evolves is similar for models and 
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Figure 9. Equivalent width of [OIII]A5007 against the mean star- 
burst age t$B for Sample I. Panel (a) shows the observed equiva- 
lent widths, whereas in (b) the dilution of W\ohi\ by the under- 
lying old stellar population is removed. Symbols as in Fig. ^. 



Figure 10. The evolution of gas excitation for Sample I galaxies, 
as measured by (a) [OIII]/H/3 and (b) [OIII]/[OII]. Reddeningand 
H/3 absorption corrections were applied. Symbols as in Fig. H. 



data. Furthermore, practically all data points are bracketed 
by the models shown! Given the already discussed caveats 
affecting both axis of this figure, it would be premature to 
use this result to draw any conclusion about, say, the IMF 
in starbursts. The point to emphasize here is that, to our 
knowledge, this is the first time that predicted and observed 
values of Wh/3 are plotted against an age axis, and it is 
gratifying to see a good agreement between data and models. 



4.2 The equivalent width of [OIII] 

The popularity of Wh/3 as an age indicator stems mostly 
from its insensitivity to nebular conditions such as density, 
temperature and metallicity. Yet, the detailed photoioniza- 
tion models for evolving starbursts by Stasihska & Leitherer 
(1996) show that the equivalent width of [OIII] is also a 
powerful chronometer of starbursts for metallicities below 
solar, as is the case for most of the galaxies studied here. 
We therefore explore the behavior of Wyoui] against our 
empirical age index Isb- 

The results for Sample I are shown in Fig. [| An anti- 
correlation is clearly present. The relation gets even stronger 
after correcting W[oiii] f° r the dilution by an underlying 
population (Fig. ^Jd) . Triangles and circles, which correspond 
to the two lower Z bins, trace rather well defined sequences 
in Fig. ^3, but the more metal rich galaxies (plotted as 
squares) present a more scattered distribution. A plausible 
explanation for this larger spread is that, as discussed by 
Stasihska & Leitherer (1996), Wyoin] ceases to be a decreas- 
ing function of age as Z approaches Zq . Fig. |i^b also shows 
that the higher Z galaxies tend to have older starbursts, an 
effect which is further discussed below. 
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Figure 11. As Fig. [To| but for non-Seyfert galaxies in Sample II. 



4.3 Gas excitation 

Another prediction of evolutionary synthesis plus photoion- 
ization calculations is that the gas excitation decreases with 
time, because the ratio of ionizing photons per gas parti- 
cle (the "ionization parameter") decreases and the ionizing 
spectrum softens as the hotter stars die (e.g., Copetti etal. 
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1986; Cid Fernandes et al. 1992; Garcia- Vargas et al. 1995; 
Stasinska & Leitherer 1996). In Figs, [l^ and |ll] we explore 
the evolution of two line ratios in order to address this issue. 

Unlike for emission line equivalent widths, particularly 
Wh/3, metallicity plays a key role in defining ratios involving 
forbidden lines because of its influence on the gas tempera- 
ture. The three Z intervals represented by different symbols 
in these figures help to disentangle the effects of evolution 
and Z. Since HII galaxies are less chemically evolved than 
Starburst nuclei, differences in metallicity should also be- 
come apparent distinguishing objects by their activity class. 
This can be readily seen in Figs. fcHJ-21. Most HII galax- 
ies (open symbols) are in the low Z bin (triangles), while 
most Starburst nuclei (filled symbols) are in the high Z bin 
(squares). Sample I contains only one Starburst nucleus of 
the 16 sources with (O/H) < 0.4 (O/H) and only 3 HII 
galaxies among the 14 objects with (O/H) > 0.6 (O/H)©. 
For Sample II, the only two Starburst nuclei are also the 
most metal rich objects. Furthermore, the four (O/H) = 
0.4-0.6 (O/H)© sources (open circles) located below the HII 
galaxy sequences in the Z-dependent Figs. [TT| a and b are 
precisely the four intermediate Starburst/HII galaxy groups 
defined by Raimann et al. (2000a) . Metallicity and activity 
class are hence practically equivalent quantities. 

Figs. ^ and [tl[show the behavior of [OIIIJ/H/3 and 
[OIII]/[OII] against tsB for Samples I and II respectively. 
In this section we discuss only results for metal poor objects 
(triangles and circles), which are mostly HII galaxies. These 
systems present clear trends of decreasing excitation for in- 
creasing Isb, in qualitative agreement with theoretical pre- 
dictions. Sources in Sample II join smoothly the sequences 
defined by sources in Sample I in all plots above, extend- 
ing it to smaller tsB- As already explained, this happens 
mainly because of its objective prism selection, which favors 
the detection of young starbursts (e.g. Stasinska & Leitherer 
1996). Since emission lines are powered solely by the most 
massive stars, they should be insensitive to the presence of 
older, non-ionizing populations, and hence we can expect 
the decrease in gas excitation to level off for Isb > 10 7 yr. 
This is consistent with the distributions of low Z objects 
in Figs, [lo] and [nj. We have also investigated other line ra- 
tios, such as [OII]/H/3 and [Nil] /Ha, both of which increase 
systematically with increasing Isb- 

These same trends, were identified and discussed by 
Stasinska etal. (2001), who used Wh/3 as an age indicator. 
This agreement is hardly surprising, since we have empiri- 
cally verified that Wh/3 and our evolutionary index Isb are 
related. In fact, as a corollary of this relation, we can auto- 
matically subscribe all trends found using Wh/3 as measure 
of evolution! We therefore need not repeat here the exten- 
sive discussions on the evolution of emission line properties 
of starbursts by Stasinska et al. (2001) and previous studies. 
The usual caveats about reddening sensitive line ratios (such 
as [OIII]/[OII]) and the effects of an absorption component 
in H/3, discussed in the references above, also apply here. 

Of course, Jsb is a new age indicator, based entirely on 
measured stellar properties. Sure enough, it too has its lim- 
itations, but these are of a completely different nature than 
the uncertainties affecting emission line age diagnostics (e.g., 
the dilution correction for Wh/3 or Wyoui})- This reassur- 
ing agreement supports the interpretation of the trends in 
[OIII]/H/3 and [OIII]/[OII] against t S B for metal poor ob- 
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Figure 12. Gas metallicity against the evolutionary index tsB- 
Symbols as in Fig. Starburst nuclei (filled symbols) are richer in 
metals than HII galaxies, and their mixed populations of 10 6 -10 8 
yr stars make them look more evolved, producing the concentra- 
tion towards to the top right. 



jects as a result of (theoretically expected) evolution of the 
gas excitation. 



4.4 Metallicity effects 

Galaxies in our (O/H) > 0.6 (O/H)© metal-rich bin, are 
heavily concentrated in the bottom-right regions of Figs. 
|l2] , corresponding to low emission line equivalent widths, 
low excitation and large age. Despite their horizontal off- 
set towards large Isb, high Z systems (plotted as squares) 
are well mixed with metal poor systems in the Wh/3 dia- 
grams (Figs. |^ and Q), in agreement with the idea that Wh/3 
is largely insensitive to Z. In the other diagrams, however, 
high Z galaxies are clearly offset along the vertical axis, 
particularly in the gas-excitation plots (Figs, [j^ and |ll[ ). 
Furthermore, while metal poor objects line up on broad but 
well defined sequences of decreasing Wiqiii], [OIIIJ/H/3 and 
[OIII]/[OII] for increasing Jsb, no clear trends appear when 
considering high Z objects by themselves. The scattered dis- 
tribution of metal rich objects in these plots is in qualitative 
agreement with the models by Stasinska etal. (2001), which 
show that Z-dependent indices such as those used in Figs. 
pl| are not good chronometers for metal rich starbursts. 

A more intriguing result is the systematic displacement 
of high Z galaxies towards large Isb- A naive interpretation 
of this offset would be that these systems represent the late 
stages of evolution of metal poor starbursts. In this scenario, 
a triangle would become a circle and then a square in Figs. fy- 
[n]. Note, however, that this would require that (O/H) in- 
creases by factors of 3-10 in less than 10 8 yr. Furthermore, 
this scenario would not be general, since there are several 
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metal poor objects with large Isb- A more appropriate read- 
ing of Figs, pj-llll is that there is a wide spread in Z for 
evolved starbursts, but there is practically no young metal 
rich system (the single exception being group G_Mrk710). 
This dichotomy is illustrated in Fig. h2l where Z is plotted 
against Isb for both samples. 

We attribute this behavior to the fact that, as already 
explained, metal rich sources are predominantly Starburst 
nuclei, whereas HII galaxies dominate the lower Z bins, as 
can be seen comparing the location of filled and open sym- 
bols in Fig. [l^. In most HII galaxies the starburst population 
is dominated by the young est ge nerations, partly due to se- 
lection effects (see Section 4.1.1) and partly due to the fact 
that these are small galaxies, where a single burst can have 
a large impact. In fact, instantaneous burst models often 
provide an acceptable description of these systems, provided 
allowance is made for the presence of an old underlying pop- 
ulation (Mas-Hesse & Kunth 1999). Starburst nuclei, on the 
other hand, present a more even distribution of stellar ages 
m the 10 6 -10 8 yr range (Langon etal. 2001), more compat- 
ible with an extended star formation episode than with a 
coeval burst. This age mixture is detected by the synthesis, 
resulting in a skew of the Isb index towards larger values. 

We therefore conclude that the trend of Z with Isb 
simply reflects the fact that metal rich objects have a more 
complex recent history of star formation than metal poor 
objects. 



5 SUMMARY 

We have investigated the evolution of emission line prop- 
erties in Starburst nuclei and HII galaxies using age diag- 
nostics based on their observed integrated stellar population 
properties. Our main results can be divided in two parts. 

In the first part of this paper, we have presented the 
results of an empirical population synthesis (EPS) analysis 
of star-forming galaxies and explored ways to condense these 
results onto simple diagrams and indices designed to assess 
the evolutionary state of a stellar population. Two useful 
tools were developed: 

(1) An evolutionary diagram: A compact description of stel- 
lar populations in terms of young (< 10 7 yr), intermediate 
age (10 s yr) and old (> 10 9 yr) components allows the evo- 
lutionary state of a galaxy to be assessed by its location on 
a [xy,Xi,xo) diagram, each axis carrying the contribution 
of stars within a given age range to the total flux. 

(2) Mean ages: Flux-weighted mean ages of both the total 
stellar population (I) and the starburst component (Isb) 
were defined . 

Both tools were tested with theoretical galaxy spec- 
tra for instantaneous bursts and continuous star formation. 
These tests showed that the evolution of stellar populations 
is adequately mapped by these empirical tools, supporting 
their application to real galaxies. Perhaps the main conclu- 
sion here is that one can achieve a good first order descrip- 
tion of the evolutionary state of a starburst using very lit- 
tle spectral information; our analysis used just 3 absorption 
lines plus 2 continuum colors in the 3600-4500 A interval. 

The EPS-analysis of two samples of starbursting galax- 
ies showed them be distributed along the direction of evolu- 



tion in the {xy ,xi ,xo) diagram. This result encouraged us 
to use our mean starburst age Isb as an empirical clock to 
gauge the evolutionary state of starbursts. 

In the second part of this study we have investigated 
correlations between the emission line properties of Star- 
burst nuclei and HII galaxies and the Isb index in order to 
test, in a completely empirical way, whether emission lines 
evolve along with the stars in starbursts. The results of this 
investigation can be summarized as follows. 

(1) We have verified that the equivalent widths of H/3 and 
[OIII] decrease for increasing ts b ■ This is in accordance with 
well known, but little tested, theoretical expectations. 

(2) The use of Wh& an d W\oiii] as age indicators is ham- 
pered by the the diluting effects of an old underlying stellar 
population unrelated to the starburst. Besides providing a 
quantitative assessment of evolution, the EPS analysis pro- 
vides a straight forward estimate of this effect. 

(3) As a whole, Starburst nuclei are found to have a more 
even distribution of stellar ages in the 10 6 -10 8 yr range than 
HII galaxies, which are often dominated by the youngest 
generations. 

(4) Three Seyfert 2 objects were also analysed, two of 
which have stellar population characteristics radically dif- 
ferent from those in starburst galaxies, as seen, for instance, 
by their location on the evolutionary diagram. These two 
sources also have dilution-corrected Wh0 values well above 
those of starbursts. The third object has characteristics sug- 
gestive of a composite starburst + Seyfert 2 system. 

(5) The gas excitation, as measured by emission line ra- 
tios, was found to decrease systematically for increasing Isb, 
also in agreement with theoretical predictions. This evolu- 
tionary sequence is only well defined for metal poor objects, 
which are mostly HII galaxies. Metal rich galaxies do not 
present clear evolutionary trends in the gas excitation in- 
dices, in qualitative agreement with photoionization models 
for evolving starbursts. 
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ABSTRACT 

In systems undergoing starbursts the evolution of the young stellar population is ex- 
pected to drive changes in the emission line properties. This evolution is usually stud- 
ied theoretically, with a combination of evolutionary synthesis models for the spectral 
energy distribution of starbursts and photoionization calculations. In this paper we 
present a more empirical approach to this issue. We apply empirical population syn- 
thesis techniques to samples of Starburst and HII galaxies in order to measure their 
evolutionary state and correlate the results with their emission line properties. A cou- 
ple of useful tools are introduced which greatly facilitate the interpretation of the 
synthesis: (1) an evolutionary diagram, whose axis are the strengths of the young, 
intermediate age and old components of the stellar population mix, and (2) the mean 
age of stars associated with the starburst, Isb- These tools are tested with grids of 
theoretical galaxy spectra and found to work very well even when only a small num- 
ber of observed properties (absorption line equivalent widths and continuum colors) 
is used in the synthesis. 

Starburst nuclei and HII galaxies are found to lie on a well defined sequence in the 
evolutionary diagram. Using the empirically defined mean starburst age in conjunction 
with emission line data we have verified that the equivalent widths of H/3 and [OIII] 
decrease for increasing tsB- The same evolutionary trend was identified for line ratios 
indicative of the gas excitation, although no clear trend was identified for metal rich 
systems. All these results are in excellent agreement with long known, but little tested, 
theoretical expectations. 

Key words: galaxies: starburst - galaxies: evolution - galaxies: stellar content - ISM: 
HII Regions 



1 INTRODUCTION 

Starburst systems in the local universe, from giant HII 
regions to Starburst nuclei, are important laboratories to 
study the evolution of massive stars and physical processes 
thought to be associated with the very early stages of galaxy 
formation. These motivations, coupled to the advances in 
modeling and observational capabilities, have led to a burst 
of activity in this field during the past decade. 

By far the most common approach to infer the physi- 
cal properties of starbursts is to compare their spectral en- 
ergy distribution (SED) with models based on evolutionary 
synthesis (eg. Mas-Hesse & Kunth 1991; Olofsson 1995; Lei- 
therer et al. 1999). This technique performs ab initio calcula- 
tions of the spectral evolution of a stellar population on the 
basis of evolutionary tracks for stars covering a wide range 
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of masses, stellar spectral libraries plus prescriptions for the 
initial mass function, star formation rate and chemical evo- 
lution. The comparison between observations and models 
may focus on particular spectral features, such as stellar 
wind lines in the UV (Robert, Leitherer & Heckman 1993), 
WR features (Cervino & Mas-Hesse 1994; Schaerer & Vacca 
1998), Balmer absorption lines (Gonzalez Delgado, Leitherer 
& Heckman 1999), lines from red super giants (Garci'a- 
Vargas et al. 1997; Mayya 1997), or on a combination of lines 
and the multiwavelength continuum (Mas-Hesse & Kunth 
1999; Lancon et al. 2001). A common difficulty faced in such 
studies is the contamination of the spectrum by an under- 
lying old stellar population, which can be significant in the 
optical-near IR range. This contamination is sometimes re- 
moved by adopting a template spectrum for the old popu- 
lation (Lancon etal. 2001), or else its effects are evaluated 
from the excess flux between models and data (Mas-Hesse 
& Kunth 1999). Other uncertainties include those associ- 
ated with differences between different sets of evolutionary 
tracks, incomplete or imperfect spectral libraries; stochas- 
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tic effects and numerical issues (Cervino etal. 2001); and 
possible effects of binary stars (Mas-Hesse & Cervino 1999). 
We refer the reader to the review by Schaerer (2001) for a 
detailed discussion. 

The massive, hot stars in young starbursts photoionize 
the surrounding gas, producing an emission line spectrum 
which can also be used to constrain the properties of the 
starburst, like its age, metallicity and star-formation rate. 
In fact, emission line diagnostics of starbursts play a central 
role in this field for the simple practical reason that emis- 
sion lines are much easier to measure than stellar features. 
A great number of studies have been devoted to develop- 
ing such diagnostics through a combination of theoretical 
SEDs from evolutionary synthesis and phototoionization cal- 
culations for the corresponding nebular conditions (Garci'a- 
Vargas, Bressan & Diaz 1995; Stasiriska & Leitherer 1996; 
Chariot & Longhetti 2001; Moy, Rocca-Volmerange & Fioc 
2001; Stasiriska, Schaerer & Leitherer 2001; Schaerer 2000). 

Among the many diagnostic tools developed, the equiv- 
alent width of H/3 (Wh@) stands out as a powerful age in- 
dicator of starbursts. The prediction, known for more than 
20 years (Dottori 1981), is that Wh& decreases as a burst 
evolves. This behavior has been so extensively confirmed by 
more elaborate calculations that Wh& is often used as a sub- 
stitute for an age axis in studies which investigate the evo- 
lution of starbursts (e.g., Stasiriska etal. 2001). Other gen- 
eral predictions of evolutionary synthesis + photoionization 
models are that the gas excitation and that the equivalent 
width of [OIII]A5007 should decrease as a starburst evolves 
(Coppeti, Dottori & Pastoriza 1986; Stasiriska & Leitherer 
1996) although such diagnostics are more critically affected 
by metallicity effects. 

These predictions, now routinely applied to infer phys- 
ical properties of starbursts, are hard to be tested directly, 
since that requires evaluating the age of a starburst without 
resorting to emission line diagnostics. In this paper we take 
a step back in time and investigate the empirical validity of 
these long known predictions by means of a simple Empir- 
ical Population Synthesis (EPS) analysis of Starburst and 
HII galaxies. This analysis allows a quantitative assessment 
of the evolutionary state of a stellar population based only 
on observed stellar features. EPS techniques have their own 
limitations (eg, Cid Fernandes etal. 2001a), but these are of 
a different nature than the uncertainties involved in evolu- 
tionary synthesis, and thus serve as an independent test of 
predictions of evolutionary synthesis models. 

Our main goals are to: 

(i) Study the stellar population properties of a large and 
varied sample of star-forming galaxies by means of an EPS 
analysis. 

(ii) Develop and test EPS-based tools to assess the evolu- 
tionary state of star-forming galaxies in a quantitative and 
easy-to-interpret manner. 

(iii) Perform a model-independent investigation of the rela- 
tion between gaseous properties and the evolutionary state 
of starbursts. 

In section 2 we present the data sets used in this study. 
Section 3 deals with points (%) and (ii) above. The EPS 
method is presented and its results for star-forming galaxies 
are discussed by means of simple empirical tools designed 
to aid the interpretation of the synthesis. We also present a 



comparative study of EPS and evolutionary synthesis meth- 
ods, which serves to test and calibrate our tools to measure 
evolution. In section 4 we address point (iii) by studying 
the evolution of emission line properties of Starburst and 
HII galaxies using an EPS-based measure of the burst age. 
Section 5 summarizes our main results. 



2 DATA 

This investigation requires optical spectra where both 
gaseous (emission lines) and stellar properties (continuum 
and absorption lines) can be discerned. Two data sets meet- 
ing this requirement were used in this work. 

The first set, which we denote "Sample I", comes from 
the studies of Storchi-Bergmann, Kinney & Challis (1995) 
and McQuade, Kinney & Calzetti (1995), extracted from the 
atlas of Kinney etal. (1993). This sample covers both large, 
luminous galaxies with star-forming nuclei ("Starburst nu- 
clei" ) and smaller, weaker systems such as HII galaxies and 
blue compact dwarves. Active galaxies were discarded, with 
the exception of NGC 6221, which is dominated by Star- 
burst activity except in X-rays (Levenson etal. 2001). We 
further limit our analysis to those objects with metallicity 
estimates by Storchi-Bergmann, Calzetti & Kinney (1994), 
which leaves a total of 41 galaxies, 18 of which are classified 
as Starburst nuclei. The spectra were collected through a 
large 10" x 20" aperture, which corresponds to a circular 
aperture of 1.3 kpc in radius at the median distance of 
the galaxies. 

Emission line fluxes and equivalent widths (W) for this 
sample were re-measured from the publicly available original 
spectra and found to be in good agreement with those ob- 
tained by Storchi-Bergmann et al. (1995) and McQuade et al. 
(1995). All spectra were corrected for Galactic extinction us- 
ing the reddening law of Cardelli, Clayton & Mathis (1989, 
with Rv = 3.1) and the E(B — V) values from Schlegel, 
Finkbeiner & Davis (1998) as listed in NED 1 . Corrections 
for internal extinction were applied based on the Ha/H/3 ra- 
tio, whose intrinsic value was taken to be 2.86, and allowing 
for underlying absorption components (see Section 4.1.1). 
For the stellar population analysis we have measured the 
W's of the Call K A3933, CN A4200 and G band A4301 
with respect to a pseudo continuum defined at selected pivot 
points, located at A = 3660, 3780, 4020 and 4510 A, fol- 
lowing the methodology outlined in Cid Fernandes, Storchi- 
Bergmann & Schmitt (1998). Our values agree very well with 
those of Storchi-Bergmann etal. (1995), but there were sig- 
nificant discrepancies between our measurements and those 
published by McQuade etal. (1995). 

"Sample II" comes from the Spectrophotometric At- 
las of HII galaxies of Terlevich etal. (1991), as analysed 
by Raimann etal. (2000a,b). Most of the individual spectra 
in this atlas do not have enough signal to measure stellar 
features, which prompted Raimann et al. to average them 
in order to increase the stellar signal. Out of 185 galaxies, 
they have defined 19 groups of similar characteristics. Each 

1 The NASA/IPAC Extragalactic Database (NED) is operated 
by the Jet Propulsion Laboratory, California Institute of Tech- 
nology, under contract with the National Aeronautics and Space 
Administration. 



© 0000 RAS, MNRAS 000, 000-000 



The evolution of stars and gas in starburst galaxies 3 



group is then treated as if it corresponded to an individ- 
ual galaxy. Three of these groups are composed of Seyfert 2 
galaxies. These were kept in our analysis only to illustrate 
their systematically different properties with respect to the 
remaining groups, 10 of which are composed of HII galax- 
ies, 2 are Starburst nuclei and 4 are classified as intermediate 
HH/Staburst systems. The typical aperture covered by these 
spectra correspond to an equivalent radius of 0.8/ifj 1 kpc. 

Raimann et al. (2000a) have measured absorption line 
Ws and continuum fluxes for these groups following the 
same methodology as above. Emission line properties were 
analysed by Raimann et al. (2000b). Line fluxes were initially 
measured after subtraction of a stellar population model 
which included internal reddening. The resulting Ha/H/3 
ratio was used to further correct for residual extinction to- 
wards the line emitting regions. We have adopted both stel- 
lar and nebular properties as published by these authors 
without further corrections. 

In summary, of all stellar and nebular properties com- 
piled for these 2 samples, the following will be used in the 
analysis below: (1) the W's of the Call K, CN and G-band 
absorption features; (2) continuum fluxes at 3600, 4020 and 
4510 A; (3) emission line fluxes and Ws of strong opti- 
cal lines: [OII]A3727, H/3, [OIII]A5007, Ha and [NII]A6584; 
(4) nebular abundances, as listed in Storchi-Bergmann et al. 
(1994) and Raimann et al. (2000b); (5) the "activity class", 
as reported in the afore mentioned papers. This last item is 
used to distinguish small systems like HII and blue compact 
galaxies from Starburst nuclei, which live on larger and more 
luminous galaxies, usually late type spirals. The latter kind 
of galaxies present a more complex mixture of stellar popu- 
lations and are richer in heavy elements than HII galaxies, 
as will become clear in the analysis that follows. 



3 EMPIRICAL POPULATION SYNTHESIS 
ANALYSIS 

3.1 The method: input and output quantities 

In order to provide a quantitative description of the stel- 
lar populations for galaxies in Samples I and II, we have 
used their absorption line Ws and continuum colors as in- 
put to the EPS algorithm developed by Cid Fernandes et al. 
(2001a). The code decomposes a spectrum onto a base of 12 
simple stellar populations of different ages and metallicities 
(Z). This base was defined by Schmidt et al. (1991) out of a 
large sample of star clusters originally observed by Bica & 
Alloin (1986a,b). The main output of the code is the pop- 
ulation vector x, whose 12 components carry the fractional 
contributions of each base element to the observed flux at a 
given normalization wavelength Ao- This vector corresponds 
to the mean solution found from a 10 8 steps likelihood- 
guided Metropolis walk through the parameter space. Since 
colors are also modeled, extinction enters as an extra param- 
eter, but this will not be directly used in our analysis. Some 
EPS studies (e.g., Bica 1988) impose that solutions follow 
well behaved paths on the age-Z plane spanned by the base, 
in order to force consistency with simple scenarios for chem- 
ical evolution. Here we follow Schmidt et al. (1991) in not 
imposing such a priori constraints in order to allow for more 
general scenarios, such as systems undergoing mergers. 



There are no major conceptual differences between this 
EPS method and that originally developed by Bica (1988) 
or its variants, which have been applied to many stellar pop- 
ulation studies (e.g., Bica, Alloin & Schmidt 1990; De Mello 
et al. 1995; Kong & Cheng 1999; Schmitt, Storchi-Bergmann 
& Cid Fernandes 1999; Raimann etal. 2000a). However, in 
this work we will explore novel ways of expressing the re- 
sults of the synthesis, which use the population vector to 
construct easy-to-interpret diagrams and indices. 

The results presented below were all obtained feeding 
the EPS code with just 5 observables: The Ws of Call K, 
CN and the G-band, plus the F3660/-F4020 and F4510/F4020 
continuum colors. The errors on these quantities were fixed 
at 0.5 A for W K and W G , 1 A for W C jv, and 0.05 for the 
colors. As discussed by Cid Fernandes etal. (2001a), the 
combination of observational errors, little input information 
and quasi-linear dependences within the base hinders accu- 
rate estimates of all 12 components of x, but reliable results 
are obtained grouping x components of same age. We have 
therefore employed age-grouping schemes in our analysis. 

The base spans 5 logarithmicaly spaced age bins: 10 6 , 
10 7 , 10 8 , 10 9 and 10 10 yr. Components with these ages are 
combined onto a reduced 5-D population vector, whose com- 
ponents are denoted by xq, x-j, x&, xq and xio respectively. 
We will also work with an even further reduced descrip- 
tion of stellar populations, in which the 10 9 and 10 10 yr 
old components are re-grouped onto xo = xq + xio, and 
the young 10 6 and 10 7 yr components are binned onto 
xy = x% + X7. Renaming the "intermediate age" 10 8 yr 
component to xi = x&, we obtain a compact, 3-D version of 
the population vector: (xy , xi , xo)- 

Normalization requires that xy + xi + xo = 1, while the 
positivity constraint implies that all components are > 0. 
Therefore, an EPS solution is confined to a triangular cut 
of a plane in the (xy , xi, xo) space, which facilitates the 
visualization of results (see Cid Fernandes et al. 2001b for 
an application of this scheme to trace the evolution of cir- 
cumnuclear starbursts in active galaxies). We note that the 
description of a galaxy spectrum in terms of x-components 
is analogous to a Principal Component Analysis (e.g., Sodre 
& Stasihska 1999), with the difference that, by construction, 
each component has a known physical meaning. 

We concentrate our analysis of the EPS results on evo- 
lutionary effects. Furthermore, we focus on the evolution of 
recent stellar generations (< 10 8 yr), associated with the 
star-forming activity in Starburst and HII galaxies. Metal- 
licity effects are discussed using the nebular oxygen abun- 
dance, which essentially reflects the metallicity of the most 
recent stellar generation (Storchi-Bergmann etal. 1994). 



3.2 EPS results and the evolutionary diagram 

The results of the EPS analysis of galaxies in Samples I 
and II are listed in Table 1 both for the (xq, x-j, x&, xq, xio) 
and the Young, Intermediate and Old descriptions. The nor- 
malization wavelength is Ao = 4020 A. A very convenient 
way to present the results of the synthesis is to project the 
(xy , xi , xo) vector onto the xy-xj plane. This is done in 
Figs, la and b for Samples I and II respectively. Dotted lines 
in these plots mark lines of constant xo, computed from the 
xy + xi + xo = 1 condition. Note that these are actually 
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EPS Results for Samples I and II 
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Q 
O 


6 




Q 
O 


35 ± 3 


55 




Q 
O 


41 




2 


6.3 ± 0.3 


G_Camll48-2020 


91 ± 4 


6 ± 4 


1 


± 


1 


1 


± 


1 


2 


± 1 


97 


± 


1 


3 


± 


1 


6.1 ± 0.1 


G.UM461 


84 ± 7 


11 ± 7 


1 


± 


1 


2 


± 


1 


3 


± 1 


95 


± 


2 


4 


± 


2 


6.1 ± 0.1 


G_Toll924-416 


59 ± 11 


29 ± 12 


5 


± 


3 


3 




2 


4 


± 2 


88 




4 


7 




3 


6.4 ± 0.3 


G_NGC1487 


56 ± 10 


29 ± 12 


8 


± 


4 


2 


± 


2 


4 


± 3 


85 


± 


5 


7 


± 


3 


6.5 ± 0.3 


G_Toll004-296 


47 ± 11 


40 ± 12 


4 


± 


3 


4 




2 


5 


± 3 


87 


± 


4 


9 




3 


6.5 ± 0.3 


G.UM488 


42 ± 11 


42 ± 12 


6 


± 


4 


A 


± 


2 


6 


± 3 


84 


± 


5 


9 


± 


3 


6.6 ± 0.3 


G_Tol0440-381 


45 ± 9 


23 ± 11 


22 


± 


6 


4 




2 


7 


± 4 


67 




6 


11 




4 


6.7 ± 0.4 


G.UM504 


35 ± 9 


21 ± 10 


25 


± 


7 


6 


± 


3 


13 


± 7 


56 


± 


6 


19 


± 


6 


6.9 ± 0.6 


G.UM71 


22 ± 8 


19 ± 9 


38 


± 


7 


7 




4 


14 


± 7 


41 




6 


21 




7 


7.2 ± 0.7 


G_NGC1510 


16 ± 7 


17 ± 9 


51 


± 


8 


4 


± 


3 


12 


± 7 


33 


± 


6 


16 


± 


7 


7.4 ± 0.6 


G_Cam0949-2126 


20 ± 8 


21 ± 9 


30 


± 


7 


13 




5 


17 


± 6 


41 




6 


29 




6 


7.1 ± 0.7 


G_Mrk711 


24 ± 10 


38 ± 12 


15 


± 


6 


8 


± 


4 


15 


± 7 


62 


± 


7 


23 


± 


6 


6.9 rt 0.7 


G.UM140 


25 ± 8 


16 ± 9 


44 


± 


7 


4 




3 


11 


± 6 


41 




6 


15 




6 


7.2 ± 0.6 


G.NGC3089 


20 ± 7 


12 ± 7 


44 


± 


8 


5 


± 


3 


19 


± 9 


32 


± 


6 


24 


± 


8 


7.3 ± 0.8 


G.Mrk710 


58 ± 11 


30 ± 12 


4 




3 


3 




2 


5 


± 3 


87 


± 


4 


8 




3 


6.4 ± 0.3 


G.UM477 


27 ± 8 


20 ± 10 


31 


± 


8 


6 


± 


3 


15 


± 7 


47 


± 


7 


21 


± 


7 


7.1 ± 0.7 


G.UM103 


6 ± 5 


10 ± 6 


45 


± 


8 


15 




7 


23 


± 8 


17 




6 


38 




8 


7.6 ± 1.0 


G.NGC4507 


3 ± 2 


6 ± 4 


4 


± 


3 


26 


± 


9 


61 


± 9 


9 


± 


4 


87 


± 


4 


7.1 ± 2.3 


G.NGC3281 


1 ± 1 


4 ± 2 


5 




3 


39 ± 11 


51 ± 11 


5 




3 


90 




4 


7.3 ± 2.9 



Table 1. Columns 2-6: Population vector in the (xe, X7, xg, xg, x\o) description. Columns 7, 4 and 8: xy, xi and xo respectively. All x 
components are in percentages of the total flux at Ao — 4020 A. Column 9: mean starburst age as defined in Section 3.4. Objects starting 
with a G_ are the spectral groups of Sample II. The three last entries in the table correspond to Seyfert 2 systems. 
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Figure 1. Results of the empirical population synthesis analysis for Samples I (a) and II (b), condensed on an evolutionary diagram. 



10 8 yr stars is plotted along the vertical axis. A third perpendicular axis carries the contribution xo of older populations (> 10 9 yr). 
Dotted lines indicate lines of constant xo, as labeled. HII galaxies are plotted as triangles and Starburst nuclei as squares. Filled symbols 
correspond to galaxies with WR features. Crosses in panel (b) indicate Seyfert 2s. Note that the Starburst nucleus G_Mrk710 and the 
HII galaxy G.Toll004.296, both of which show WR features, overlap at (ly.ij) = (0.87,0.04) in panel (b). 



straight lines in the xy-xi plane, which appear curved be- 
cause of the logarithmic axis. 

Galaxies from both samples define a smooth sequence 
from large xy to large xi, with not much spread in xo, par- 
ticularly for Sample II. The larger spread seen in Fig. la 
is partly due to the fact that the data for Sample I was 
collected through apertures typically 2.6 times larger than 
for Sample II, and thus sample a more heterogeneous mix 
of stellar populations. This interpretation is supported by 
the fact that NGC 1510, which appears in both samples, 
looks somewhat younger in Sample II (see Table 1). Another 
source of scatter in Fig. la stems from the large number of 
Starburst nuclei in Sample I. These systems, represented 
by squares in both panels, live on galaxies with a signifi- 
cant old stellar component, whose effect is to drag points 
towards the bottom- left of the plot. For instance, the point 
at (xY,xi,xo) = (0.28,0.06,0.66) in Fig. la is NGC 5996, 
whose spectrum reveals a weak starburst immersed in an 
old population (McQuade etal. 1995; Kennicutt 1992). HII 
galaxies, on the other hand, are closer to "pure starbursts" . 
In fact, they were once thought to be truly young galaxies 
undergoing their first star-formation episode (Searle & Sar- 
gent 1972). Only recently it has become clear that they too 
contain old stars (Telles & Terlevich 1997; Schulte-Ladbeck 
& Crone 1998; Raimann etal. 2000a). This explains why 
Sample II, which is essentially an HII galaxy sample, exhibits 
a more well defined sequence in Fig. 1, with all non-AGN 
sources bracketed by the xo = and 30% contours. 



The three deviant crosses spoiling the Sample II se- 
quence in Fig. lb are the Seyfert 2 groups, with their pre- 
dominantly old stellar populations (Raimann etal. 2000a). 
G-UM103 has a significant "post-starburst" component, 
reminiscent of more evolved starburst + Seyfert 2 composite 
systems, while the other two groups occupy a region charac- 
teristically populated by LINERs and non-composite Seyfert 
2's (Cid Fernandes etal. 2001b). 

Since the location of a galaxy in Fig. 1 reflects the evo- 
lutionary state of its stellar population, we interpret the dis- 
tribution of objects in this diagram as an evolutionary se- 
quence, with the mean stellar age running counter-clockwise. 

There are several reasons to interpret Fig. 1 as an evo- 
lutionary sequence. First, metal absorption lines become 
deeper and galaxy colors become progressively redder as one 
moves from large xy to large xi along the sequence. In fact, 
the sequence defined by Sample II follows very closely the 
blue to red (young to old) spectral sequence in Figure 1 of 
Raimann etal. (2000a). Second, all Sample II galaxies in 
which WR features have been detected (those marked by 
filled symbols in Fig. lb) are located in the large xy re- 
gion of the diagram, consistent with the young burst ages 
(a few Myr) implied by the mere presence of WR stars. 
Filled symbols in Fig. la mark galaxies from Sample I 
which are listed in the WR-galaxy catalog maintained by 
D. Schaerer (webast.ast.obs-mip.fr/people/scharer). Their 
more even distribution, as compared to Sample II, is due 
to the old population and aperture effects discussed above 
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and nicely illustrated by Meurer (2000). Whereas the data 
analyzed here pertains to kpc-scales, spectra used to classify 
a starburst as a WR-galaxy are usually obtained through 
much narrower slits centered on the brightest cluster, thus 
favoring the detection of young systems. Processing such 
spectra through our EPS-machinery would surely move the 
filled points in Fig. la towards younger ages. (Conversely, 
one would expect that narrow slit observations of galaxies 
represented by empty symbols in the bottom right of Fig. la, 
such as Mrk 357 and UGCA 410, have a good chance of re- 
vealing WR features.) Finally, galaxies located in the large 
xi zone in the top-left (such as NGC 1800 in Sample I and 
group G-NGC3089 in Sample II) have spectra typical of 
a "post-starburst" population, with pronounced high order 
Balmer absorption lines typical of A stars (Gonzalez Del- 
gado etal. 1999). It is important to remark that neither the 
presence of WR features nor Balmer absorption lines were 
used in the EPS analysis, and yet the EPS results are com- 
patible with the information apported by these observables. 

3.3 EPS analysis of theoretical galaxy spectra 

A straight-forward theoretical reason to interpret Fig. 1 as 
an evolutionary sequence is that, schematicaly, a simple 
(i.e., coeval) stellar population moves on this diagram from 
(xy,xi,xo) = (1,0,0) at age t = to (0,1,0) after some 
10 8 yr and then to (0, 0, 1) for ages > 10 9 yr. 

In order to follow this evolutionary path more closely 
we have carried out an EPS analysis of theoretical galaxy 
spectra from GISSEL96, the evolutionary synthesis code of 
Bruzual & Chariot (1993). The theoretical spectra were pro- 
cessed in exactly the same way as the real spectra of Samples 
I and II. Instantaneous burst and continuous star-formation 
models were computed for various ages between t = and 15 
Gyr, a Salpeter IMF between 0.1 and 125 M© and Z = Z Q . 
GISSEL96 uses stellar tracks from the Padova group and 
offers a choice of spectral libraries. We have chosen the one 
which uses the atlas of Jacoby, Hunter & Christian (1984) for 
the optical range. The absorption features necessary for our 
EPS are clearly defined with this library. Spectral resolution 
was in fact the reason we have chosen GISSEL96 over the 
Starburst99 code of Leitherer etal. (1999), which is more 
taylored to study young stellar populations but currently 
works with an optical library too coarse for EPS analysis. 



3.3.1 Burst models 

The results for an instantaneous burst are shown as a solid 
line in the (iy,i;, xo) evolutionary diagram of Fig. 2a. La- 
bels next to selected points indicate the model age in Myr. 
As expected, evolution proceeds from xy to xi to xo, such 
that a position on this diagram can be associated with an 
age. Fig. 3a provides an alternative representation of this 
diagram, in which all 3 components are explicitly plotted in 
a face-on projection of the plane containing the (xy , xi, xo) 
vector. The idea for this projection was borrowed from simi- 
lar plots by Pelat (1997, 1998) and Moultaka & Pelat (2001). 

In principle one would expect all evolution prior to ages 
of ~ 10 9 yr to progress along the xo = contour, whereas 
in practice the GISSEL96 models oscillate between xo = 10 
and 15% for t < 2 x 10 8 yr. Similarly, the value of xj starts 



to increase before t = 10 7 yr, when all stars should still be- 
long to the xy age bin. These deviations occur due to the 
limited number of observables used in the synthesis and be- 
cause these contain observational errors which broaden the 
likelihood-function of x in a non-trivial way (Cid Fernandes 
et al. 2001a). As a result, some of the true xy proportion al- 
ways spills over onto xj and xo, and so on. Overall, however, 
these figures show an excellent correspondance between em- 
pirical and evolutionary populations synthesis calculations. 

In Fig. 4c we show the behavior of the 5 age components 
xg-x io as a function of the age of the GISSEL96 models. 
The population vector evolves smoothly with age, except 
for the small kink a little short of 10 7 yr due to the sudden 
appearance of red supergiants (Chariot & Bruzual 1991). 
As expected, x% peaks around t = 10 6 yr, X7 peaks around 
t = 10 7 yr and so on, but note how X7 is less well defined 
than any other component. The EPS decomposition tends 
to represent a ~ 10 7 yr burst as a combination of xq and xg 
instead of a strong X7. Also, a non-negligible fraction of the 
10 6 yr component spills from x% onto X7 for t ^ 10 6 yr. A 
similar effect occurs with xq and xio for t <; 10 yr. 

Such imprccisions in the mapping between evolutionary 
and empirical population synthesis are largely suppressed 
in the coarser, but more robust, (xy , xi, xo) description, 
as shown in Fig. 4b. The figure also shows the evolution of 
xsb = x^+X7+xs, which we hereafter treat as the "starburst 
component" , representing the past ~ 10 8 yr of the history 
of star formation in a galaxy. This is a more reasonable 
definition for our purposes than using only the youngest, 
ionizing population (xq), since single burst models are not 
adequate to describe kpc scale regions such as those sampled 
by the observations of Samples I and II. Instead, the inner 
kpc of star-forming galaxies, particularly Starburst nuclei, 
contains a collection of many individual associations plus a 
field population with a spread in age. The detailed studies 
by Langon etal. (2001) and Tremonti etal. (2001) illustrate 
this point (see also Calzetti 1997; Legrand etal. 2001). Such 
systems are frequently better represented by models with 
multiple bursts or continuous star formation over ~ 10 8 yr 
(Meurer 2000; Meurer 1995; Coziol, Barth & Demers 1995; 
Coziol, Doyon & Demers 2001). 

3.3.2 Burst plus an underlying old population 

Real galaxies have a mixture of stellar populations of differ- 
ent ages, and galaxies in Samples I and II are no exception. 
The ongoing star-formation which makes them classifiable 
as starburst systems is observed atop an old (> 10 9 yr) stel- 
lar substrate formed in the earlier history of the galaxy. In 
our young, intermediate and old description, the effect of 
this underlying population is to dilute the values of xy and 
xi, which represent the recent history of star-formation. As 
a result, an instantaneous burst occurring on top of an old 
background does not follow the evolutionary sequence traced 
by the solid line in Fig. 2a. 

Two quantities suffice to examine these diluting effects: 
the fraction xoid(0) of the total Lsurst + Lou luminosity at 
the start of the burst (t = 0) which is due to old stars, and 
the function l(t) = LBurst(t)/LBurst(0), which describes the 
luminosity evolution of the burst in units of its initial lumi- 
nosity. Naturally, all these quantities refer to the same wave- 
length, Ao. With these definitions, and considering that Lou 
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Figure 2. Evolution of GISSEL96 models on the (xy ,xj ,xo) EPS diagram, obtained processing the theoretical galaxies as the data. 
Dotted lines indicate lines of constant xo- (a) The solid line shows results for an instantaneous burst. Numbers indicate the model age in 
Myr (note that evolution proceeds counterclockwise). Dashed lines correspond to instantaneous bursts on top of an old population which 
at t — accounts for xoidi®) — 10% ( t0 P) an d 50% (bottom) of the flux at Aq = 4020 A. (b) Models with continuous star formation. 
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does not evolve significantly on the time-scales of interest 
(< 10 s yr), it is easy to show that 



6 7 8 9 

log model age [yr] 

Figure 4. EPS parameters against model age for instantaneous 
burst models, (a) Mean ages, (b) population vector x in the re- 
duced Young, intermediate and old description, and (c) in the 
5 component description. The Isb curve in (a) is stopped at a 
few times 10 s yr, after which it loses meaning for instantaneous 
bursts. 



xoid{t) 



XBurst(0)l(t) + XQld(0) 



(1) 



where XBurst{0) = 1 — »ow(0). We can now look at the evo- 
lution of x presented above for a pure burst as corresponding 
just to the LBurst(t) component, which at time t accounts 
for only XBurst(t) = 1 — xoid{t) of the total luminosity of 
the system. This allows us to, with the help of equation 1, 
re- normalize the evolution of (xy , xi , xo) to this new scale 
for any desired value of the contrast parameter xoid(O)- 

Results for xoid(0) = 10 and 50% are shown as dashed 
lines in Figs. 2a and 3a. As the burst fades, the evolution- 
ary sequence bends over towards large xo quicker for larger 
xoid(0), i.e., for smaller initial ratios of burst to underly- 
ing old population power. Despite this effect, evolution still 
proceeds in an orderly counterclockwise fashion. 



3.3.3 Continuous star-formation models 

Figs. 2b and 3b show the evolution of [xy , xi, Xo) for GIS- 
SEL96 models with continuous star-formation. Since in this 
regime there are always young stars, at any given age > the 
system looks younger than an instantaneous burst model. At 
t = 10 s yr, for instance, (xy,xi,xo) ~ (0.6,0.3,0.1) for the 
continuous star-formation models and (0.3,0.6,0.1) for an 
instantaneous burst. For t > 10 9 yr, x converges to a re- 
gion around (0.3, 0.4, 0.3), instead of plunging towards large 
xo as for instantaneous bursts. Evolutionary sequences for 
models with time-decaying star-formation rates (say, expo- 
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Figure 3. Evolution of GISSEL96 models in the (xy ,xj ,xo) EPS-space seen from a face on projection of the xy + xj + xq = 1 plane, 
(a) Instantaneous burst models. The top curve corresponds to a pure burst, while the middle and bottom curves correspond to bursts 
atop a pre-existing old population with xoidi®) = 10 an d 50% respectively, (b) Continuous star formation models. 



nentially or an "extended burst" step function), would define 
curves intermediary between those traced in Figs. 2a and b. 

Since the luminosity increases without bounds for con- 
tinuous star formation, any underlying old population is 
quickly outshone by the new stars. We therefore do not 
present dilution curves such as those computed for an in- 
stantaneous burst, since, except for the earliest ages, they 



are practically identical to the undiluted curves in Figs. 2b 
and 3b. The evolution of the population vector for contin- 
uous star-formation models is shown in Figs. 5b and c. As 
expected, the curves are smoother, and each component lives 
longer than for an instantaneous burst (Fig. 4). 
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6X6 + 7X7 + 8X8 



(4) 
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log model age [yr] 

Figure 5. As Fig. 4 for models with continuous star formation. 



3.4 Mean stellar age 

The experiments above demonstrate that the evolutionary 
state of a starburst can be assessed by its location on the 
(xy,xi,xo) diagram. In order to translate this location 
into a number which quantifies the "evolutionary state" one 
could use, for instance, the angle 6 = tan _1 x//xy, which 
increases as a burst evolves. Alternatively, we may use x 
to compute the mean age t of the stellar population. Since 
stellar populations evolve in a non-linear way, it makes more 
sense to define t from the mean log t among the populations 
represented by the base: 



logi(Ao) = ^ Xj(Xo) log tj 



(2) 



The dependence of x on the normalization wavelength 
is explicitly written in this equation to emphasize that this 
definition of? is A-dependent. This happens because the Xi's 
are flux fractions at Ao, so ?(Ao) is a flux-weighted mean age. 
A A-dependent age makes observational sense for the sim- 
ple reason that young stars are bluer than old stars, which 
makes t an increasing function of A. Though the value of t 
depends on Ao, the evolutionary sequence traced by this in- 
dex independs on the choice of normalization, and so it can 
be used to rank populations on different evolutionary states. 

For our 5-ages base (i< = 10 s , 10 7 , 10 8 , 10 9 and 10 10 
yr), ?(Ao) becomes (in yr) 



log t(Xo) = 6X6 + 7x7 + 8xs + 9xg + lOxio 



(3) 



Since we are primarily interested in quantifying the evo- 
lutionary stage of populations associated with the most re- 
cent star-formation in starburst systems, it is interesting to 
consider a definition of t which removes the diluting effects 
of an underlying old stellar population. This can be done 
re- normalizing X6 + X7 + xs to 1, which yields the following 
definition for the mean starburst age (also in yr): 



Xe + X7 + Xg 

Note that by construction 6 < log? < 10 and 6 < logtsB < 
8. 

The solid lines in Fig. 4a compare our t and Isb EPS- 
based age indices for Ao = 4020 A with the corresponding 
age of the GISSEL96 models for an instantaneous burst. 
Despite some minor oscillations, these two indices bear a ~ 
one-to-one relation with the theoretical age. Fig. 5a presents 
these same indices but for the continuous star-formation 
models. As expected, t and ts b evolve more slowly than for 
an instantaneous burst, but they still increase steadily with 
the model age. Since t and !sb are entirely obtained from a 
few easily measurable quantities, this result encourages their 
use as empirical clocks for stellar populations. 

As it is clear from its very definition, due to the coarse 
age-resolution of the base, our mean age index Isb is not 
meant to be used as a fine-graded chronometer of starbursts. 
Yet, the above experiments with theoretical spectra clearly 
show that Isb provides a useful way to rank galaxies ac- 
cording to the age of the dominant population among the 
multiple generations of stars formed in the recent history 
of star-formation. This definition is particularly well suited 
to describe spectra which are integrated over large regions 
and hence average over many such generations, Although 
all galaxies discussed here contain populations younger than 
10 7 yr, which power their emission line spectrum, this ongo- 
ing star-formation may be less intense than in the recent past 
(10 7 -10 8 yr), such that the young generations live among an 
older, non-ionizing starburst population. In this case, one 
expects to find significant X7 and xs components, and thus 
tsB > 10 7 yr. Conversely, if the current star-formation is 
more vigorous than in the past, mean ages of less than 10 7 
yr are expected. It is in this context of starbursts extended 
over a period of up to ~ 10 s yr that we envisage Isb and 
the (x.y , xi , xo) diagram as useful tracers of evolution. 

In principle, a base with a finer age resolution, including 
elements intermediate between X6, X7 and xs, could yield 
a more detailed description of the evolution of starbursts. 
In practice, however, these elements would be well approxi- 
mated by linear combinations of the existing base elements 
unless new observables were introduced in the synthesis pro- 
cess. For this reason, we opted to perform our EPS analy- 
sis with the base and observables described in §3.1, whose 
pros and cons have already been fully exploited in our pre- 
vious investigations (Cid Fernandes et al. 2001a,b; Schmitt 
et al. 1999). Furthermore, as we shall soon see, this relatively 
coarse description is well suited to our present purposes. 



4 THE EVOLUTION OF EMISSION LINE 
PROPERTIES 

Emission lines in star-forming galaxies are umbilicaly linked 
to their young stellar population, whose massive, hot stars 
photoionize the surrounding gas. Also, in non-instantaneous 
starbursts the continuum carries a large contribution of stars 
borne before the current generation of ionizing stars, thus 
affecting emission line equivalent widths. In this section, we 
combine the tools to measure the evolution of starbursts de- 
veloped in Section 3 with the emission line data compiled in 
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Figure 6. Equivalent width of H/3 against the EPS-defined mean 
starburst age tsB for Sample I. Panel (a) shows the observed 
equivalent widths (W^|). In (b) the contribution of old stars to 
the continuum under H/3 is removed dividing by £sb(4861). 

In panel (c) Wh@ is further corrected for the presence of an ab- 
sorption component. Triangles, circles and squares correspond to 
gas metallicity ranges of (O/H) < 0.4, 0.4-0.6 and > 0.6 (O/H) 
respectively. Open and filled symbols are used to distinguish HII 
galaxies from Starburst nuclei. The P r values are the probabili- 
ties of no correlation in a Spearman test, small values indicating 
significant correlations. 



Section 2 to investigate whether the emission line properties 
do indeed evolve along with the burst. 



Figure 7. As Fig. 6 but for Sample II. Crosses indicate Seyfert 
2's, which were excluded from the correlation analysis. 



are determined from completely independent measurements, 
which only highlights the significance of this result. 

The different symbols in Figs. 6, 7 and all subsequent 
plots represent three gas metallicity ranges, triangles, cir- 
cles and squares corresponding to (O/H) < 0.4, 0.4 to 0.6 
and > 0.6 (O/H) respectively. Open and filled symbols are 
used to distinguish HII galaxies from Starburst nuclei. We 
postpone a discussion of the effects of Z and activity class to 
Sections 4.3 and 4.4, which explore emission line properties 
more directly affected by these variables. 



4.1 The equivalent width of H/3 

As a burst ages, its most massive stars are the first to die, 
resulting in a steady decline of the ionizing photon flux 
and hence on the luminosity of recombination lines such 
as H/3. The stellar continuum Ch/3 underneath H/3 also de- 
creases, but more slowly than L H p{t), since it carries a sig- 
nificant contribution from longer-lived, non-ionizing, lower 
mass stars. As a result Wap = Lh^/Cmp decreases as the 
burst evolves, as first discussed by Dottori (1981) and con- 
firmed by evolutionary synthesis calculations. Wap is there- 
fore an age indicator, and it is frequently used as such in 
studies of star-forming systems (e.g, Stashiska et al. 2001). 
However, an empirical confirmation of the Wupit) predic- 
tion requires an independent measure of the starburst age. 

In Figs. 6a and 7a we carry out this test with galaxies 
from Samples I and II respectively, using our EPS-based 
index tsa as a clock for the starburst. The anti-correlation 
between Wgp and Isb is evident for both samples, thus 
confirming that Wufi decreases with time. The probability 
P r of no correlation in a Spearman's rank test is just 7x 10 -5 
for Sample I and 10~ 4 for Sample II, indicating a very high 
statistical significance. We emphasize that Wh@ and !sb 



4-1.1 Corrections to Wnp 

The values of Wnp in Figs. 6a and 7a are the raw measure- 
ments (Wff b ^). At least two corrections have to be consid- 
ered, both of which increase Wh@- 

(1) The continuum under H/3 carries a contribution from 
an old stellar population which dilutes Wh/3 with respect to 
the value it would have in a pure starburst. Our EPS analysis 
provides a natural way of correcting for this effect, which is a 
major source of concern in studies which use Wh/3 as an age- 
indicator (Stasihska et al. 2001 and references therein). In or- 
der to isolate the contribution of the starburst to Ch/3 it suf- 
fices to multiply it by «sb(4861), the starburst component 
defined in Section 3.3.1 but renormalized to A = 4861 A. The 
corrected Wh/3 is thus simply IFJ^ /xsb(4861). Figs. 6b and 
7b show the dilution-corrected evolution of Wh@- 

The effects of this correction are largest for two of the 
Seyfert 2 groups in Sample II, which move significantly 
above the sequence in Fig. 7b because of their bulge dom- 
inated optical continuum (large xo)- The effect is not so 
large for G_UM103, which, as already discussed, resembles 
a starburst + Seyfert 2 composite. The fact that the ob- 
served values of Wh@ in Seyfert 2's are in general smaller 
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than those in starburst systems is purely due to this dilution 
effect. As explained by Cid Fernandes etal. (2001b), Seyfert 
2's ought to have intrinsically larger Wh@ than starbursts, 
as we obtain with our EPS-based dilution correction. 

Among the star-forming galaxies in Samples I and II, 
this correction typically increases Wh/3 by ~ 50% (in good 
agreement with the corrections inferred by Mas-Hesse & 
Kunth 1999 on similar objects), but it reachs more than 
a factor of 2 in some cases. The correction is somewhat 
smaller for Sample II, partly because it contains intrinsically 
younger systems and partly because of its smaller apertures, 
which reduces "contamination" by an extended old stellar 
population. The dilution correction improves the Wh0-Jsb 
correlation for Sample I and degrades the one for Sample 
II, while for the combined sample the statistical significance 
remains unchanged at P r = 5 x 10~ 10 . 

(2) A second correction to be considered is that due 
to the presence of an absorption component hidden under- 
neath the H/3 emission. This component is present in our 
spectral base with strengths of up to WJj'g = 8 A, achieved 
for 10 s yr populations. We have used the population vector 
x obtained in the synthesis to compute the expected value 
of W^, typically 3-5 A. Adding to Wgp yields a cor- 

rected emission Wh/3- The combined dilution and absorption 
corrected values of Wh/3 are shown in Figs. 6c and 7c. 

The absorption correction is only significant for galaxies 
with weak H/3 emission, such as those bellow W$ ~ 10 A 
at the bottom right of Fig. 6a. For these systems one has to 
look at the absorption-corrected values of Wh/3 as uncertain 
by as much as a factor of 2. The correction is negligible for 
most galaxies in Sample II, which, due to its objective prism 
selection criterium, contains more strong lined objects than 
Sample I. Indeed, the mean W$f) is 55 A for Sample II but 
just 22 A for Sample I. This is also why, as a whole, Sample 
II contains a higher proportion of young starbursts. 

The absorption correction degrades the Wh/3-Jsb cor- 
relation for Sample I slightly. For the combined Sample I + 
II data the P r value increases from 5 x 10~ 10 to 3 x 10~ 8 , 
which is still significant at the 5-a level. We thus see that 
these '1 st order corrections' introduce very little scatter, and 
do not alter our conclusion that Wh/3 does indeed evolve 
along with the stars that make up a starburst. One can also 
look at this result the other way around, and conclude that 
the fact that the Whp versus Jsb diagram behaves as ex- 
pected proves the usefulness of our EPS-based evolutionary 
index Jsb, with the advantage that it is immune to the di- 
lution and absorption effects which plague Wh@ and other 
emission-line based age indicators. 

We note in passing that these corrections alone are 
enough to bring the values of Wh/3 within the range spanned 
by evolutionary synthesis calculations such as those by Lei- 
therer etal. (1999), whereas, as it has long been known, 
the raw observed values fall bellow such predictions (Bres- 
solin, Kennicutt & Garnett 1999 and references therein). 
Differential extinction, with line emitting regions being more 
reddened than the stellar continuum (Calzetti, Kinney & 
Storchi-Bergmann 1994) and leakage of ionizing photons out 
of the HII regions associated with the starburst are further 
examples of processes that act in the sense of reducing Whp- 
We thus concur with Raimann et al. (2000b) and Stasiriska 
et al. (2001) in that the apparent mismatch between theo- 
retical and observed values of Wh/3 bears no fundamental 
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Figure 8. Observed and predicted evolution of Wjj/3- Data points 
are as in Figs. 6c and 7c. Lines correspond to theoretical Wga{t) 
curves with t converted onto the EPS-based Jsb scale. Solid line: 
GISSEL96 models with a Salpeter IMF between 0.1 and 125 
Mq. Dashed lines: Starburst99 models with a Salpeter IMF with 
M Mpp = 100 Mq (long dashes) and 30 M© (short dashes). 

physical significance. In fact, the combined effects of differ- 
ential extinction, leakage and the uncertainties in the dilu- 
tion and absorption corrections are probably responsible for 
most of the vertical scatter in Figs. 6c and 7c. 



4-1-2 Comparison with models 

A quantitative comparison of the predicted evolution of 
Wh& with that detected in Figs. 6 and 7 demands process- 
ing model spectra through the same EPS machinery used 
to analyse the data. This is necessary to translate model 
ages (t) onto our Jsb age scale. We have used the evolution 
of the rate of ionizing photons N(H ) and the continuum 
under H/3 predicted by GISSEL96 to compute Wnpit) for 
the same models synthesized in Section 3.3, for which the 
t ~^Jsb conversion is shown in Figs. 4a and 5a. 

The result for the continuous star formation models is 
shown as a solid line in Fig. 8, over-plotted onto the data 
points from Samples I and II (Figs. 6c and 7c). The dashed 
curves in this plot are the Starburst99 predictions (Leitherer 
et al. 1999) for solar metallicity, constant star formation 
models with a Salpeter IMF up to M upp = 100 and 30 M . 
These latter curves are drawn with the t —¥ Jsb conversion 
obtained for GISSEL96, since, as already explained, it is not 
currently possible to do an EPS analysis with Starburst99 
due to its poor spectral resolution optical libraries. Instan- 
taneous burst models (not shown for clarity) follow roughly 
the same curves up to Jsb ~ 10 6 ' 8 yr and then plunge verti- 
cally, signaling the end of the ionizing phase of the cluster. 

The rate at which Whp evolves is similar for models and 
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Figure 9. Equivalent width of [OIII]A5007 against the mean star- 
burst age Isb for Sample I. Panel (a) shows the observed equiva- 
lent widths, whereas in (b) the dilution of Wyoiii] by the under- 
lying old stellar population is removed. Symbols as in Fig. 6. 



Figure 10. The evolution of gas excitation for Sample I galaxies, 
as measured by (a) [OIII]/H/3 and (b) [OIII]/[OII]. Reddening and 
H/3 absorption corrections were applied. Symbols as in Fig. 6. 



data. Furthermore, practically all data points are bracketed 
by the models shown! Given the already discussed caveats 
affecting both axis of this figure, it would be premature to 
use this result to draw any conclusion about, say, the IMF 
in starbursts. The point to emphasize here is that, to our 
knowledge, this is the first time that predicted and observed 
values of Wup are plotted against an age axis, and it is 
gratifying to see a good agreement between data and models. 



4.2 The equivalent width of [OIII] 

The popularity of Wh/3 as an age indicator stems mostly 
from its insensitivity to nebular conditions such as density, 
temperature and metallicity. Yet, the detailed photoioniza- 
tion models for evolving starbursts by Stasiriska & Leitherer 
(1996) show that the equivalent width of [OIII] is also a 
powerful chronometer of starbursts for metallicities below 
solar, as is the case for most of the galaxies studied here. 
We therefore explore the behavior of W[oju] against our 
empirical age index tsB- 

The results for Sample I are shown in Fig. 9. An anti- 
correlation is clearly present. The relation gets even stronger 
after correcting Wtoui] for the dilution by an underlying 
population (Fig. 9b). Triangles and circles, which correspond 
to the two lower Z bins, trace rather well defined sequences 
in Fig. 9b, but the more metal rich galaxies (plotted as 
squares) present a more scattered distribution. A plausible 
explanation for this larger spread is that, as discussed by 
Stasihska & Leitherer (1996), Wtoui] ceases to be a decreas- 
ing function of age as Z approaches Zq . Fig. 9b also shows 
that the higher Z galaxies tend to have older starbursts, an 
effect which is further discussed below. 
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Figure 11. As Fig. 10 but for non-Seyfert galaxies in Sample II. 



4.3 Gas excitation 

Another prediction of evolutionary synthesis plus photoion- 
ization calculations is that the gas excitation decreases with 
time, because the ratio of ionizing photons per gas parti- 
cle (the "ionization parameter" ) decreases and the ionizing 
spectrum softens as the hotter stars die (e.g., Copetti et al. 
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1986; Cid Fernandes etal. 1992; Garda-Vargas et al. 1995; 
Stasinska & Leitherer 1996). In Figs. 10 and 11 we explore 
the evolution of two line ratios in order to address this issue. 

Unlike for emission line equivalent widths, particularly 
Whp, metallicity plays a key role in defining ratios involving 
forbidden lines because of its influence on the gas tempera- 
ture. The three Z intervals represented by different symbols 
in these figures help to disentangle the effects of evolution 
and Z. Since HII galaxies are less chemically evolved than 
Starburst nuclei, differences in metallicity should also be- 
come apparent distinguishing objects by their activity class. 
This can be readily seen in Figs. 6-12. Most HII galax- 
ies (open symbols) are in the low Z bin (triangles), while 
most Starburst nuclei (filled symbols) are in the high Z bin 
(squares). Sample I contains only one Starburst nucleus of 
the 16 sources with (O/H) < 0.4 (O/H) and only 3 HII 
galaxies among the 14 objects with (O/H) > 0.6 (O/H)©. 
For Sample II, the only two Starburst nuclei are also the 
most metal rich objects. Furthermore, the four (O/H) = 
0.4—0.6 (O/H)© sources (open circles) located below the HII 
galaxy sequences in the Z-dependent Figs. 11a and b are 
precisely the four intermediate Starburst/HII galaxy groups 
defined by Raimann etal. (2000a). Metallicity and activity 
class are hence practically equivalent quantities. 

Figs. 10 and 11 show the behavior of [OIIIJ/H/3 and 
[OIII]/[OII] against Isb for Samples I and II respectively. 
In this section we discuss only results for metal poor objects 
(triangles and circles), which are mostly HII galaxies. These 
systems present clear trends of decreasing excitation for in- 
creasing Isb , in qualitative agreement with theoretical pre- 
dictions. Sources in Sample II join smoothly the sequences 
defined by sources in Sample I in all plots above, extend- 
ing it to smaller Isb- As already explained, this happens 
mainly because of its objective prism selection, which favors 
the detection of young starbursts (e.g. Stasinska & Leitherer 
1996). Since emission lines are powered solely by the most 
massive stars, they should be insensitive to the presence of 
older, non-ionizing populations, and hence we can expect 
the decrease in gas excitation to level off for Isb > 10 7 yr. 
This is consistent with the distributions of low Z objects 
in Figs. 10 and 11. We have also investigated other line ra- 
tios, such as [OIIJ/H/3 and [NII]/Ha, both of which increase 
systematically with increasing Isb- 

These same trends, were identified and discussed by 
Stasinska etal. (2001), who used Wh/3 as an age indicator. 
This agreement is hardly surprising, since we have empiri- 
cally verified that Wnp and our evolutionary index tsB are 
related. In fact, as a corollary of this relation, we can auto- 
matically subscribe all trends found using Wh/3 as measure 
of evolution! We therefore need not repeat here the exten- 
sive discussions on the evolution of emission line properties 
of starbursts by Stasinska et al. (2001) and previous studies. 
The usual caveats about reddening sensitive line ratios (such 
as [OIII]/[OII]) and the effects of an absorption component 
in H/3, discussed in the references above, also apply here. 

Of course, Isb is a new age indicator, based entirely on 
measured stellar properties. Sure enough, it too has its lim- 
itations, but these are of a completely different nature than 
the uncertainties affecting emission line age diagnostics (e.g., 
the dilution correction for Wh& or W[ ni})- This reassur- 
ing agreement supports the interpretation of the trends in 
[OIII]/H/3 and [OIII]/[OII] against t SB for metal poor ob- 
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Figure 12. Gas metallicity against the evolutionary index Isb- 
Symbols as in Fig. 6. Starburst nuclei (filled symbols) are richer in 
metals than HII galaxies, and their mixed populations of 10 6 -10 8 
yr stars make them look more evolved, producing the concentra- 
tion towards to the top right. 

jects as a result of (theoretically expected) evolution of the 
gas excitation. 

4.4 Metallicity effects 

Galaxies in our (O/H) > 0.6 (O/H)© metal-rich bin, are 
heavily concentrated in the bottom-right regions of Figs. 6- 
12, corresponding to low emission line equivalent widths, 
low excitation and large age. Despite their horizontal off- 
set towards large Isb, high Z systems (plotted as squares) 
are well mixed with metal poor systems in the Wjjp dia- 
grams (Figs. 6 and 7), in agreement with the idea that Wh/3 
is largely insensitive to Z. In the other diagrams, however, 
high Z galaxies are clearly offset along the vertical axis, 
particularly in the gas- excitation plots (Figs. 10 and 11). 
Furthermore, while metal poor objects line up on broad but 
well defined sequences of decreasing W[ m], [OIIIJ/H/3 and 
[OIII]/[OII] for increasing tsB, no clear trends appear when 
considering high Z objects by themselves. The scattered dis- 
tribution of metal rich objects in these plots is in qualitative 
agreement with the models by Stasinska et al. (2001), which 
show that Z-dependent indices such as those used in Figs. 6- 
11 are not good chronometers for metal rich starbursts. 

A more intriguing result is the systematic displacement 
of high Z galaxies towards large tsB- A naive interpretation 
of this offset would be that these systems represent the late 
stages of evolution of metal poor starbursts. In this scenario, 
a triangle would become a circle and then a square in Figs. 6- 
11. Note, however, that this would require that (O/H) in- 
creases by factors of 3-10 in less than 10 s yr. Furthermore, 
this scenario would not be general, since there are several 
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metal poor objects with large Isb- A more appropriate read- 
ing of Figs. 6-11 is that there is a wide spread in Z for 
evolved starbursts, but there is practically no young metal 
rich system (the single exception being group G_Mrk710). 
This dichotomy is illustrated in Fig. 12, where Z is plotted 
against Isb for both samples. 

We attribute this behavior to the fact that, as already 
explained, metal rich sources are predominantly Starburst 
nuclei, whereas HII galaxies dominate the lower Z bins, as 
can be seen comparing the location of filled and open sym- 
bols in Fig. 12. In most HII galaxies the starburst population 
is dominated by the youngest generations, partly due to se- 
lection effects (see Section 4.1.1) and partly due to the fact 
that these are small galaxies, where a single burst can have 
a large impact. In fact, instantaneous burst models often 
provide an acceptable description of these systems, provided 
allowance is made for the presence of an old underlying pop- 
ulation (Mas-Hesse & Kunth 1999). Starburst nuclei, on the 
other hand, present a more even distribution of stellar ages 
in the 10 6 -10 8 yr range (Lancon etal. 2001), more compat- 
ible with an extended star formation episode than with a 
coeval burst. This age mixture is detected by the synthesis, 
resulting in a skew of the Isb index towards larger values. 

We therefore conclude that the trend of Z with Isb 
simply reflects the fact that metal rich objects have a more 
complex recent history of star formation than metal poor 
objects. 



5 SUMMARY 

We have investigated the evolution of emission line prop- 
erties in Starburst nuclei and HII galaxies using age diag- 
nostics based on their observed integrated stellar population 
properties. Our main results can be divided in two parts. 

In the first part of this paper, we have presented the 
results of an empirical population synthesis (EPS) analysis 
of star-forming galaxies and explored ways to condense these 
results onto simple diagrams and indices designed to assess 
the evolutionary state of a stellar population. Two useful 
tools were developed: 

(1) An evolutionary diagram: A compact description of stel- 
lar populations in terms of young (< 10 7 yr), intermediate 
age (10 8 yr) and old (> 10 9 yr) components allows the evo- 
lutionary state of a galaxy to be assessed by its location on 
a (xy , xi , xo) diagram, each axis carrying the contribution 
of stars within a given age range to the total flux. 

(2) Mean ages: Flux-weighted mean ages of both the total 
stellar population (t) and the starburst component (Isb) 
were defined . 

Both tools were tested with theoretical galaxy spec- 
tra for instantaneous bursts and continuous star formation. 
These tests showed that the evolution of stellar populations 
is adequately mapped by these empirical tools, supporting 
their application to real galaxies. Perhaps the main conclu- 
sion here is that one can achieve a good first order descrip- 
tion of the evolutionary state of a starburst using very lit- 
tle spectral information; our analysis used just 3 absorption 
lines plus 2 continuum colors in the 3600-4500 A interval. 

The EPS-analysis of two samples of starbursting galax- 
ies showed them be distributed along the direction of evolu- 



tion in the (xy , xi, xo) diagram. This result encouraged us 
to use our mean starburst age Isb as an empirical clock to 
gauge the evolutionary state of starbursts. 

In the second part of this study we have investigated 
correlations between the emission line properties of Star- 
burst nuclei and HII galaxies and the Isb index in order to 
test, in a completely empirical way, whether emission lines 
evolve along with the stars in starbursts. The results of this 
investigation can be summarized as follows. 

(1) We have verified that the equivalent widths of H/3 and 
[OIII] decrease for increasing Is b ■ This is in accordance with 
well known, but little tested, theoretical expectations. 

(2) The use of Wh@ and W[om] as a g e indicators is ham- 
pered by the the diluting effects of an old underlying stellar 
population unrelated to the starburst. Besides providing a 
quantitative assessment of evolution, the EPS analysis pro- 
vides a straight forward estimate of this effect. 

(3) As a whole, Starburst nuclei are found to have a more 
even distribution of stellar ages in the 10 6 -10 8 yr range than 
HII galaxies, which are often dominated by the youngest 
generations. 

(4) Three Seyfert 2 objects were also analysed, two of 
which have stellar population characteristics radically dif- 
ferent from those in starburst galaxies, as seen, for instance, 
by their location on the evolutionary diagram. These two 
sources also have dilution-corrected Wh/3 values well above 
those of starbursts. The third object has characteristics sug- 
gestive of a composite starburst + Seyfert 2 system. 

(5) The gas excitation, as measured by emission line ra- 
tios, was found to decrease systematically for increasing Isb, 
also in agreement with theoretical predictions. This evolu- 
tionary sequence is only well defined for metal poor objects, 
which are mostly HII galaxies. Metal rich galaxies do not 
present clear evolutionary trends in the gas excitation in- 
dices, in qualitative agreement with photoionization models 
for evolving starbursts. 
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